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1.  Introduction 


The  computer  code  SDA  to  be  described  in  this  report  solves 
the  staircase  integer  linear  prograaaing  problem,  (§>?)  given  by: 


maximize 


r 

ctxt 


subject  to:  A,  xt 


IjXj  ♦  S ^  t-2/3/(«a/I 

i  xf  i  Uj/  t=l,2,...,T 

x£  integer,  t=l,2,...,T. 


t— 1, 2, •••/!• 


This  forauiation  arises  in  aultipiant  production  allocation 
probleas,  aultisector  econoaic  planning  aodels,  and  aultitlae  period 
production  and  Inventory  problems. 

The  solution  method  used  by  the  prograa  SDA  relies  upon 
decomposition  of  the  problem  (SP)  into  smaller  subproblems  (S£(),  each 
of  which  can  then  be  efficiently  solved  by  an  LP~based 


br anch-and-bound  routine.^Each  subproblea  (S^ )  is  specified  by  a 


vector  of  costs  c{,  the  "diagonal"  submatrix  Af,  the  "of f diagonal" 
submatrix  the  right-hand  side  bg,  and  upper  and  lower  bounds  Ut 

and  Lg  on  the  variables  x^.  The  actual  formulation  is  as  follows: 


<v 


maximize 
subjec  ,  to: 


ct*t 

*t*t  *  b«  ■  B*1  V. 

Lt  *  *t  *  ut 


The  algorithm  proceeds  by  finding  a  solution  to  a  subproblea 
(say  )  and  moving  forward  to  the  next  subproblea  (St^()  with  a  new 
right-hand  side  determined  by  the  solution  to  S^.  When  a  set  of 
solutions  to  all  the  subproblems  is  found,  this  solution  will  be 
feasible  for  the  original  problem.  Prices  are  calculated  to  help 
guide  the  search  for  subproblea  solutions  toward  an  optimal  value, 
bounds  on  the  aaximua  objective  value  obtainable  given  the  current 
solutions  to  a  subset  of  the  subproblems  are  utilized  to  speed 
fathoming.  See  Pollenz  C1980]  for  a  detailed  discussion  of  the 
staircase  decomposition  algorithm. 

The  branch-and-bound  search  procedure,  which  is  used  to  solve 
the  subproblems  obtained  via  decomposition  of  (SP),  is  a  slightly 
modified  version  of  the  computer  program  BB  written  by  Gary  Kochman  as 
part  of  his  dissertation  at  Stanford  University.  This  computer  code 
solves  pure  integer  programming  problems  with  general  upper  and  lower 
bounds  on  the  variables  using  a  branch-and-bound  technique  similar  to 
that  presented  by  Dakin  C19653.  Tomlin's  (19713  improved  penalties 
are  employed  to  guide  the  choice  of  branching  variable  at  unfathomed 
nodes,  with  the  branch  being  taken  in  the  direction  opposite  to  the 
maximum  penalty.  Nodes  are  removed  from  the  branch-and-bound  list 
according  to  a  last-in-f irst-out  (LIFO)  strategy,  and  reoptimization 
after  a  branch  is  accomplished  by  the  dual  simplex  method.  See 
Kochman  C19761  for  further  details  about  this  program. 

The  linear  programming  portions  of  the  code  BB  were  developed 
by  John  Tomlin  of  the  Systems  Optimization  Laboratory  at  Stanford 
University  and  adapted  by  Kochman  to  efficiently  deal  with  simple 
upper  and  lower  bounds  on  the  variables.  The  most  important  features 
of  LPM-1  are:  storage  of  the  basis  inverse  in  product  form  (see 
Orchard-Hays  C19683);  LU  decomposition  of  the  basis  invsrse  (ssa 
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Benichou  et  al.  C19773);  and  storage  of  all  data/  Including  the 
inverse/  in  coapressed  fora  (i.e.,  zeroes  are  not  saved). 

In  order  to  aaintain  compatibility  with  BB/  the  coaputer  code 
SDA  was  written  in  FORTRAN.  Testing  and  evaluation  of  this  prograa 
was  carried  out  on  the  SCORE  DEC  Systea  20  coaputer  at  Stanford 
University.  Due  to  the  restrictive  nature  of  the  FORTRAN-10  coapiler 
at  SCORE/  feu  of  the  features  of  SDA  would  be  unacceptable  to  another 
system.  The  sole  exception  is  the  usage  of  the  systea  clock  routine 
IHPTIM/  uhich  is  employed  to  record  total  CPU  tiae  and  the  percentage 
of  execution  time  devoted  to  various  portions  of  the  algorithm.  Such 
information  is  valuable/  but  if  no  systea  clock  is  available/  the 
appropriate  sections  of  the  prograa  can  be  eliminated  without 
affecting  the  whole. 

The  main  prograa  a:.d  input  requirements  for  SOA  are  discussed 
in  the  next  section.  The  current  restrictions  on  problem  size  are 
given  in  section  3.  The  various  subroutines  which  comprise  the  bulk 
of  the  program  SDA  are  outlined  in  section  4.  The  output  generated  by 
SDA  is  described  in  section  5/  and  saaple  input  and  output  are  given 
at  the  end  of  the  appendix/  following  the  listing  of  the  prograa. 


The  main  prograa  calls  the  input  subroutine/  a  subroutine 
which  generates  a  few  initial  bounds/  and  then  iteratively  calls  the 
branch-and-bound  program  to  solve  the  appropriate  subproblems.  After 
the  optiaal  solution  has  been  discovered  (and  optiaality  verified)/ 
control  passes  to  the  next  section  of  the  main  prograa/  which  computes 
some  timing  information  and  calls  the  output  subroutine. 


The  first  input  card  aust  contain  values  for  the  parameters 


IFPROB,  NP,  IQBJ,  INVFRQ,  ITRFRQ,  and  INITBD,  in  (414,  15,  HO) 
format.  These  variables  have  the  following  significance: 


IFPROB  : 

NP 

ioaj 

INVFKQ 

ITRFRQ 

IMITBD 

read  in 


problem  identification  number  (must  be  nonzero). 

number  of  time  periods  (equivalently,  number  of  subproblems ). 

row  number  of  objective  row.  Default  is  l.  (Currently,  IOBJ 
must  be  1  for  subroutine  BOUNDR  to  operate  correctly;  this 
can  be  corrected.) 

frequency  with  which  basis  inversion  is  carried  out  in  the 
linear  programming  portions  of  the  code.  If  INVFRQ  =  k,  basis 
reinversion  will  occur  after  every  k  simplex  iterations. 
INVFRQ  should  not  be  set  greater  than  the  number  of  rows  in 
the  smallest  suhproblem. 

upper  limit  on  total  number  of  simplex  iterations. 
Computations  are  terminated  if  the  total  number  of  simplex 
iterations  exceeds  ITRFRQ.  If  ITRFRQ  =  0,  the  default  value 
of  999999  is  used. 

initial  lower  bound  estimate  for  maximum  objective  value.  It 
is  used  in  fathoming  tests  until  the  first  feasible  solution 
with  objective  value  greater  than  INITBO  is  found. 

Following  the  initial  input  card,  data  for  each  subproblee  is 
sequentially.  The  data  must  be  input  in  the  following  order: 


1.  A  leading  card  with  'SUB  nnnn'  in  columns  5-12,  where  nnnn 


is  replaced  by  the  appropriate  subproblaa  nuabor.  Subprobloa  nuabars 
aust  ba  sequent i alr  beginning  with  1. 

2.  Nonzero  entries  of  the  of  {diagonal  subaatrlx  (except  for 
subprobloa  1#  for  uhich  no  offdiagonal  aatrix  is  needed). 

3.  Rou  naae  and  type  for  each  constraint. 

4.  Nonzero  entries  of  the  diagonal  aatrix  A^. 

5.  The  right-hand  side  bt. 

6.  Lower  bounds  Lt#  followed  by  upper  bounds  U^#  in  foraat  15F5.0. 

The  type  of  constraint  is  denoted  by  a  single  letter  preceding 
the  rou  naae.  'E*  stands  for  an  equality#  *L*  and  'G*  identify  less 
than  or  equal  to  and  greater  than  or  equal  to  inequalities 
respectively#  and  aarks  the  objective  rou.  The  foraat  for  this 
rou  identification  input  is  (IX#  XI#  2X#  A8). 

For  the  reading  of  all  coefficients#  the  following  inforaation 
oust  be  specified:  coluan  naae#  rou  naae#  and  entry  value.  The  foraat 
used  for  this  is  <4X#  AS#  2X#  AS#  2X#  F12.4#  3X#  AS#  2X,  F12.4).  The 
pattern  for  rou  naae  and  value  is  repeated  so  that  two  entries  for  the 
sane  coluan  may  be  input  on  one  line*  The  rou  naaes  aust  aatch  a  naae 
read  in  earlier  (input  step  ff3).  Zero  coefficients  say  be  oaitted 
coapletely.  This  foraat  is  the  saae  as  used  for  LPM-1#  and  is 
consistent  uith  standard  MPS  foraat. 

The  only  restriction  on  input  data  is  that  the  cost 
coefficients  are  assuaed  to  be  integral.  Furtheraore#  these 
coefficients  aust  be  input  with  the  opposite  sign  (i.e.#  if  you  wish 
to  aaxiaize  cx#  you  oust  input  -c  for  the  objective  rou).  This  is  an 
unfortunate  result  of  the  fact  that  LPM-1  is  a  ainiaization  routine. 

After  all  data  has  been  read  in#  an  END  card  signifies  the  end 
of  the  data  file. 

The  input  for  a  saaple  run  of  SDA  is  Included  in  the  appendix# 


following  the  prograa  listing. 


3. 


U&£.£X-2J24 


For  current  dimensioning#  the  following  restrictions  apply  to 
use  of  the  SDA: 


1.  The  nuaber  of  subprobieas  (NP)  must  be  £  10. 

2.  The  total  nuaber  of  rows,  excluding  the  objective  row,  must  not 
exceed  60. 

3.  The  total  number  of  coluans  (including  slacks)  must  be  £  120. 

4.  The  total  number  of  nonzero  eleaents  in  the  constraint  matrix 

(excluding  the  eleaents  of  A, )  must  not  exceed  1000.  (This 
restriction  is  necessary  only  for  subroutine  BOUNDR.) 

5.  The  total  nuaoer  of  nonzero  eleaents  in  the  first  diagonal 
subaatrix#  A, #  must  not  exceed  1000. 

6.  No  right-hand  side  bt(i)  aay  exceed  1000  (due  to  default  upper 
bounds  on  the  slack  variables). 

7.  The  number  of  nodes  in  the  branch-and-bound  tree  must  never 

exceed  500.  (For  CO/ 13  problems#  the  nuaber  of  nodes  will 

never  exceed  the  total  number  of  columns#  so  this  restriction 
will  not  be  a  factor.) 


BLOCK  DATA  (from  LPM-1):  sets  initial  values  for  many  global  prograa 
constants#  Including  aaxiaua  problem  dimensions  and  ainiaua 


tolerances* 


SUBROUTINE  INPUT  (IFPROB/INITBD)  (fro*  LPM-1):  reads  In  all  problea 
data,  checks  that  problea  diaensions  do  not  exceed  current 
specif ications. 

Paraaeters: 

IFPROB  -  nonzero  problea  identification  nuaber  (output) 

INITBD  -  initial  louer  bound  estiaate  of  aaxiaal  objective 
value  (output  parameter) 

SUBROUTINE  INPSTQ:  stores  all  data  relevant  to  subproblea  whose  nuaber 
is  stored  in  variable  NS.  This  subroutine  is  called  after  reading  the 
data  initially  (from  INPUT)  and  before  each  forward  step  of  the 
algorithm. 

Entry  points: 

STORE  ~  entry  point  froa  subroutine  TESTX.  After  an  inteqral 
solution  is  found  to  subproblea  NS*  the  current  state 
of  the  corresponding  search  is  saved  by  a  call  to 
STORE. 

SUBROUTINE  RESTOR  (MNFLAG):  restores  all  data  from  subproblea  NS  in 
preparation  for  either  a  forward  or  a  backtracking  step.  For  a 
forward  step,  the  LP  relaxation  of  this  subproblea,  with  new 
right-hand  side  determined  by  a  call  to  FIXRHS,  is  solved. 

Paraaeters: 

MNFLAG  -  Input  paraaeter  set  at  0  if  backtracking,  1  if  taking 
a  forward  step. 


SUBROUTINE  FIXRHS:  commutes  new  right-hand  side  for  subproblea  NS, 
given  the  (newly  established)  setting  of  variables  of  the  previous 


subproblem 


SUBROUTINE  BOUNDR:  calculates  two  LP-based  bounds  on  the  maximum 
objective  values  for  soae  aggregations  of  the  subproblems.  These 
bounds  are  used  for  fathoming  in  subroutines  TESTX/  BKTRAK/  and 
PENLTS.  Also/  a  vector  of  prices  for  use  in  guiding  the 
branch-and-bound  search  (see  subroutine  PENLTS)  is  computed. 

SUBROUTINE  UPDATX:  updates  right-hand  side  by  taking  into 
consideration  variables  which  are  nonbasic  at  their  upper  or  lower 
(nonzero)  bounds. 

SUBROUTINE  FTRAN  (IPAP)  (from  LPM-1):  perforas  forward  transformation 
of  aatrix  coluan  (stored  in  vector  V)  by  basis  inverse. 

Parameters: 

IPAR  -  input  parameter  indicating  which  eta-vectors  are  used 
to  update  the  aatrix  column. 

SUBROUTINE  BTRAN  (from  LPM-1):  performs  backward  transformation  on 
column  stored  in  vector  ¥. 

SUBROUTINE  FORMC  (from  LPM-1):  forms  objective  function  vector  for 
Phase  1  of  primal  simplex  method. 

SUBROUTINE  PRICE  (from  LPM-1):  prices  out  nonbasic  columns  for  prlaal 
simplex  method  and  chooses  pivot  column  (stored  in  variable  JCQLP). 
Also  checks  for  dual  feasibility. 

SUBROUTINE  CHUZR  (from  LPM-1):  chooses  pivot  row  for  primal  simplex 
method  using  min  ratio  test;  stores  pivot  row  in  variable  IROWP. 
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SUBROUTINE  MRETA  (from  LPM-1):  forms  a  new  eta-vector  for  the  product 
for®  of  the  inverse. 

SUBROUTINE  SHIFTR  ( IOLD,  INEM)  (from  LPM-1):  rearranges  data  storage. 
Parameters: 

IOLD/  INEW  -  input  parameters  Indexing  storage  locations.  Move 
from  location  designated  by  IOLD  to  that 
designated  by  INEW. 

SUBROUTINE  INVERT  (from  LPM-1):  Creates  basis  inverse  by  LU 
decomposition. 

SUBROUTINE  UNPACK  (IV)  (froa  LPM-1):  expands  compressed  matrix  columns 
by  inserting  zeroes  appropriately. 

Parameters: 

IV  -  input  parameter  indexing  the  matrix  column  to  be 
expanded. 

SUBROUTINE  SHFTE  (from  LPM-1):  Subroutine  for  INVERT. 

SUBROUTINE  'JPBETA  (from  LPM-1):  updates  right-hand  sides  following  a 
primal  or  dual  simplex  pivot. 

SUBROUTINE  NORMAL  (ITSINV)  (from  LPM-1):  directs  execution  of  primal 
simplex  method. 

Parameters: 

ITSINV  -  counts  number  of  simplex  iterations  since  last  basis 
Inversion  (for  comparison  with  1NVFRQ). 
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SUBROUTINE  BANDB  ( INITfaD)  (fro*  8B):  master  program  for 
branch-and-bound  Integer  programming  routine*  It  also  serves  as 
•aster  proqram  for  reoptimization  via  the  revised  dual  siapiex  method 
after  a  forward  branch. 

Parameters: 

INITBD  -  initial  lover  bound  estimate  on  maximum  objective 
value. 

Entry  points: 

BnENTR  -  reentry  point  from  main  program  to  apply  branch-and- 
bound  search  to  any  subproblem  after  the  first. 

SUBROUTINE  DCHUZR  (from  BB):  Selects  pivot  row  for  dual  simplex  method 
and  stores  it  in  variable  IROWP.  Also  checks  for  dual  optimality. 

SUBROUTINE  DCHUZC  (from  BB):  selects  pivot  column  for  dual  simplex 
method  and  stores  it  in  variable  JCOLP.  Also  checks  for  dual 
feasibility. 

SUBROUTINE  TESTX  (from  BB):  tests  LP-optimal  solution  at  current  node 
for  integrality  and  for  fathoming  (fathoming  tests  and  branching 
strategy  modified  substantially  from  those  of  BB).  Any  new  complete 
solution  found  is  immediately  printed  out  and  saved  in  the  array 
INCUMB)  any  integral  solution  to  a  subproblem  (other  than  the  last) 
causes  a  call  to  entry  STORE  in  preparation  for  a  forward  step. 
Variable  MSTAT  flags  the  result  of  testing  and  is  checked  in 
subroutine  BANDB. 

SUBROUTINE  PENLTS  (from  BB):  computes  Tomlin's  improved  up-  and  down- 
penalties  and  the  Goaory  penalty  at  each  node.  Checks  for  forced 
branches  on  both  basic  and  non-basic  variables.  If  fathoming  does  not 


occur/  a  branch  variable  is  chosen  in  accordance  with  the  combination 
of  Tomlin's  penalties  and  the  prices  computed  in  BOUNDR/  and 
subroutine  BRANCH  is  called*  (Substantial  modifications  have  been  made 
from  the  version  of  this  subroutine  given  in  BB.) 

SUBROUTINE  BRANCH  (from  BB):  Performs  necessary  bookkeeping  for 
branching  on  variable  indexed  by  ICQL.  Increments  list  of  stored 
nodes/  revises  bounds  on  branching  variable/  and  saves  opposite  branch 
direction  and  a  bound  on  the  maximum  objective  value  on  that  opposite 
branch. 

SUBROUTINE  BKTRAK  (from  BB):  backtracks  to  a  promising  (unfathomed) 
node  from  the  list  of  stored  nodes.  Employs  last- in-first-out 
selection  rule.  If  backtracking  brings  us  back  to  the  previous 
subproblem/  the  appropriate  data  and  status  of  the  search  of  that 
subproblem  are  restored  by  a  call  to  subroutine  RESTOR. 

SUBROUTINE  MRAPUP  (from  BB):  Outputs  final  solution  information.  (See 
output  from  sample  run  at  the  end  of  the  appendix.) 

For  further  information  on  subroutines  derived  from  LPM-1  see 
Tomlin  C 19753. 

The  output  produced  by  this  proqram  falls  onto  3  sections. 

The  initial  output  contains  the  problem  identification  number  and 
dimensions/  followed  by  the  prices  and  bounds  computed  by  subroutine 
BOUNOR.  At  this  point  the  iterative  portion  of  the  algorithm  is 


begun.  Each  tine  an  improved  feasible  solution  is  discovered,  the 
tine  of  discovery,  total  nuaber  of  branches  taken,  and  new  naxiaua 
objective  value  (called  INCVAL)  are  printed  out  by  subroutine  TESTX. 
After  termination,  soae  information  regarding  the  total  tiae  taken  and 
tiae  spent  at  certain  tasks  is  output.  If  termination  occurs 
normally,  the  final  section  of  output  contains  the  optimum  solution 
and  its  objective  value.  The  solution  is  printed  in  2015  format  in 
the  following  order:  objective  value,  slack  variables,  and  integer 
variables  for  subproblem  1,  then  subproblea  2,  ...,  and  finally  for 
subproblem  NP.  If  computation  is  halted  because  the  limit  on 
iterations  (ITRFRQ)  has  been  exceeded,  then  the  best  solution  found  so 
far  ana  the  total  number  of  simplex  iterations  used  are  printed. 
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Appendix: 

Prograa  Listing*  Saaple  Input  and  Output* 
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****  snplt  Input  **** 


1  4 

1 

5  9999 

0 

SUB 

1 

ROMS 

N  OBJ 

L  ROM 

101 

L  ROM 

102 

L  ROM 

103 

L  ROM 

104 

L  ROM 

105 

COLUMNS 

COL 

101 

OBJ 

-10.00 

ROM 

101 

-10.00 

COL 

101 

ROM 

102 

-6.00 

ROM 

103 

18.00 

COL 

101 

ROM 

104 

7.00 

ROM 

105 

9.00 

COL 

102 

OBJ 

-19.00 

ROM 

101 

-2.  00 

COL 

102 

ROM 

102 

-9.00 

ROM 

103 

15.00 

COL 

102 

ROM 

104 

15.00 

ROM 

105 

0.00 

COL 

103 

OBJ 

-7.00 

ROM 

101 

4.00 

COL 

103 

ROM 

102 

9.00 

ROM 

103 

0.00 

COL 

103 

ROM 

104 

-9.00 

ROM 

105 

-8.00 

COL 

104 

OBJ 

-12.00 

ROM 

101 

9.00 

COL 

104 

ROM 

102 

15.00 

ROM 

103 

-10.00 

COL 

104 

ROM 

104 

4.00 

ROM 

105 

-8.00 

COL 

105 

OBJ 

-6.00 

ROM 

101 

5.00 

COL 

105 

ROM 

102 

15.00 

ROM 

103 

12.00 

COL 

105 

ROM 

104 

0.00 

ROM 

105 

0.00 

RHS 

RHS1 

ROM 

101 

39.00 

ROM 

102 

36.00 

RHSl 

ROM 

103 

25.00 

ROM 

104 

20.00 

RHS1 

ROM 

105 

35.00 

BOUNDS 

0. 

0. 

0. 

0.  0. 

1. 

1. 

1. 

IV  1. 

SUB 

2 

ROMS 

N  OBJ 

L  ROM 

201 

L  ROM 

202 

L  ROM 

203 

L  ROM 

204 

L  ROM 

205 

OPFDIAGQNJkL  COLUMNS 

COL 

101 

ROM 

201 

1.00 

ROM 

202 

12.00 

COL 

101 

ROM 

203 

17.00 

ROM 

204 

0.00 

COL 

101 

RON 

205 

3.00 

COL 

102 

ROM 

201 

-2.00 

ROM 

202 

0.00 

COL 

102 

ROM 

203 

16.00 

ROM 

204 

17.00 

COL 

102 

ROM 

205 

13.00 

COL 

103 

ROM 

201 

-8.00 

ROM 

202 

12.00 

COL 

103 

ROM 

203 

0.00 

ROM 

204 

-1.00 

COL 

103 

ROM 

205 

5.00 

COL 

104 

ROM 

201 

17.00 

ROM 

202 

7.00 

COL 

104 

ROM 

203 

-2.00 

ROM 

204 

-1.00 

COL 

104 

ROM 

205 

3.00 

COL 

105 

ROM 

201 

-4.00 

ROM 

202 

7.00 

COL 

105 

ROM 

203 

0.00 

ROM 

204 

0.00 

COL 

105 

ROM 

205 

-8.00 

COLUMNS 

COL 

201 

OBJ 

-4.00 

ROW 

201 

7.00 

COL 

201 

ROW 

202 

3.00 

ROW 

203 

-1.00 

COL 

201 

ROW 

204 

0.00 

ROW 

205 

8.00 

COL 

202 

OBJ 

-12.00 

ROW 

201 

6.00 

COL 

202 

ROW 

202 

0.00 

ROW 

203 

-9.00 

COL 

202 

ROW 

204 

16.00 

ROW 

205 

17.00 

COL 

203 

OBJ 

-4.00 

ROW 

201 

19.  00 

COL 

203 

ROW 

202 

-6.00 

ROW 

203 

-10.00 

COL 

203 

ROW 

204 

14.00 

ROW 

205 

20.00 

COL 

204 

OBJ 

-15.00 

ROW 

201 

12.00 

COL 

204 

ROW 

202 

15.00 

ROW 

203 

15.00 

COL 

204 

ROW 

204 

0.00 

ROW 

205 

0.00 

COL 

205 

OBJ 

-12.00 

ROW 

201 

-2.00 

COL 

205 

ROW 

202 

0.00 

ROW 

203 

2.00 

COL 

205 

ROW 

204 

8.00 

ROW 

205 

10.00 

RHS1 

ROW 

201 

31.00 

ROW 

202 

29.00 

RHS1 

ROW 

203 

26.00 

ROW 

204 

23.00 

RHS1 

ROW 

205 

23.00 

BOUNDS 

0.  0*  0.  Oa 
1.  1.  1.  la 

SUB  3 

ROMS 
N  OBJ 
L  ROW  301 
L  ROW  302 
L  ROW  303 
L  ROW  304 
L  ROW  305 
OFPDI AGONAL  COLUMNS 


COL 

201 

ROW 

301 

9.00 

ROW 

302 

COL 

201 

ROW 

303 

-9.00 

ROW 

304 

COL 

201 

ROW 

305 

6.00 

COL 

202 

ROW 

301 

0.00 

ROW 

302 

COL 

202 

ROW 

303 

8.00 

ROW 

304 

COL 

202 

ROW 

305 

11.00 

COL 

203 

ROW 

301 

18.00 

ROW 

302 

COL 

203 

ROW 

303 

-3.00 

ROW 

304 

COL 

203 

ROW 

305 

-6.00 

COL 

204 

ROW 

301 

-7.00 

ROW 

302 

COL 

204 

ROW 

303 

-10.00 

ROW 

304 

COL 

204 

ROW 

305 

-10.00 

COL 

205 

ROW 

301 

14.00 

ROW 

302 

COL 

205 

ROW 

303 

-2.00 

ROW 

304 

COL 

205 

ROW 

305 

0.00 

COLUMNS 

COL 

301 

OBJ 

-2.00 

ROW 

301 

COL 

301 

ROW 

302 

0.00 

ROW 

303 

COL 

301 

ROW 

304 

-6.00 

ROW 

305 

COL 

302 

OBJ 

-12.00 

ROW 

301 

COL 

302 

ROW 

302 

-6.00 

ROW 

303 

COL 

302 

ROW 

304 

-2.00 

RON 

305 

COL 

303 

OBJ 

-18.00 

ROW 

301 

COL 

303 

ROW 

302 

0.00 

ROW 

303 

COL 

303 

ROW 

304 

-4.00 

ROW 

305 

COL 

304 

OBJ 

-2.00 

ROW 

301 

COL 

304 

ROW 

302 

-6.00 

ROW 

303 

COL 

304 

ROW 

304 

10.00 

ROW 

305 

COL 

305 

OBJ 

-12.00 

ROW 

301 

17.00 

8.00 

15.00 

20.00 

-5.00 

14.00 

-1.00 

4.00 

0.00 
-4.  00 


9.00 

8.00 

0.00 
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COL 

305 

RON 

302 

5.00 

RON 

303 

2.00 

COL 

305 

RON 

304 

-4.00 

RON 

305 

-10.00 

RHS 

RHS 1 

RON 

301 

22.00 

RON 

302 

38.00 

RHS1 

ROW 

303 

21.00 

RON 

304 

23.00 

RHS1 

RON 

305 

21.00 

BOUNDS 

0. 

0. 

0. 

0.  0. 

1. 

1. 

1. 

1.  1. 

SUB 

4 

RONS 

N  OBJ 

L  RON 

401 

L  ROW 

402 

L  RON 

403 

L  RON 

404 

L  RON 

405 

OFFDIAGQNAL  COLUMNS 

COL 

301 

RON 

401 

1.00 

RON 

402 

-2.  00 

COL 

301 

RON 

403 

-6.00 

ROW 

404 

14.00 

COL 

301 

RON 

405 

4.00 

COL 

302 

RON 

401 

1.00 

RON 

402 

4.  00 

COL 

302 

RON 

403 

15.00 

RCN 

404 

19.00 

COL 

302 

RON 

405 

10.00 

COL 

303 

RON 

401 

-4.00 

ROW 

402 

11.00 

COL 

303 

RON 

403 

18.00 

ROW 

404 

-3.00 

COL 

303 

ROW 

405 

17.00 

COL 

304 

ROW 

401 

5.00 

RON 

402 

6.00 

COL 

304 

RON 

403 

9.00 

RON 

404 

16.00 

COL 

304 

ROW 

405 

14.00 

COL 

305 

RON 

401 

4.00 

RON 

402 

20.00 

COL 

305 

RON 

403 

10.00 

ROW 

404 

3.00 

COL 

305 

RON 

405 

16.00 

COLUMNS 

COL 

401 

OBJ 

-1.00 

RON 

401 

7.00 

COL 

401 

RON 

402 

17.00 

ROW 

403 

3.00 

COL 

401 

RON 

404 

12.00 

ROW 

405 

16.00 

COL 

402 

OBJ 

-16.00 

ROW 

401 

12.  00 

COL 

402 

RON 

402 

14.00 

ROW 

403 

11.00 

COL 

402 

ROW 

404 

12.00 

ROW 

405 

-10.00 

COL 

403 

OBJ 

-12.00 

ROW 

401 

17.00 

COL 

403 

RON 

402 

6.00 

ROW 

403 

4.00 

COL 

403 

RON 

404 

-8.00 

ROW 

405 

12.00 

COL 

404 

OBJ 

-14.00 

ROW 

401 

5.00 

COL 

404 

RON 

402 

10.00 

RON 

403 

11.00 

COL 

404 

RON 

404 

-4.00 

ROW 

405 

18.00 

COL 

405 

OBJ 

-10.00 

ROW 

401 

16.00 

COL 

405 

RON 

402 

17.00 

RON 

403 

6.00 

COL 

405 

RON 

404 

-6.00 

RON 

405 

3.00 

RHS 

RHS1 

RON 

401 

38.00 

ROW 

402 

29.00 

RHSi 

RON 

403 

29.00 

ROW 

404 

36.00 

RHS1 

RON 

405 

27.00 

BOUNDS 

0. 

0. 

0. 

0.  0. 

1. 

1. 

1. 

1.  1. 

tor 
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****  Saaple  Output  **** 
PROBLEM  1 

20  ROWS  20  COLUMNS  4  PERIODS 

PRICE  OF  OPFDIAGONAL  COLUMNS  IN  SUB.  4 
-2.29  4.57  12.57  6.86  22.86 

PRICE  OF  OFFDIAGONAL  COLUMNS  IN  SUB.  3 

0.00  0.00  0.00  0.00  0.00 

PRICE  OF  OFFDIAGONAL  COLUMNS  IN  SUB.  2 
2.12  9.18  3.53  2.12  -5.65 

MAXIMUM  CUMULATIVE  OBJECTIVE  VALUES 
FOR  LAMBDA  =  .5:  196  140  89 

FOR  LAMBDA  =  0:  90  49  42 


INTERMEDIATE  SOLUTIONS  FOUND 


TIME 

= 

0.90 

seconds; 

BRANCHES  = 

21; 

INCVAL  = 

80 

TIME 

s 

0.97 

SECONDS; 

BRANCHES  = 

28; 

INCVAL  * 

84 

TIME 

0.99 

SECONDS; 

BRANCHES  = 

31; 

INCVAL  = 

90 

time 

s 

1.32 

SECONDS; 

BRANCHES  = 

76; 

INCVAL  = 

92 

TIME 

s 

1.84 

SECONDS; 

BRANCHES  = 

145; 

INCVAL  = 

105 

TOTAL  SOLUTION  TIME  =  2.12  SECONDS 

TIME  FOR  INPUT  =  0.44  SECONDS 

TIME  FOR  LP  SOLUTIONS  =  0.31  SECONDS 

TIME  FOR  BOUNDING  &  PRICING  =  0.32  SECONDS 

TIME  FOR  BOOKKEEPING  OPERATIONS  =  0.33  SECONDS 

NUMBER  OF  FORWARD  STEPS  TAKEN  =  80 

TOTAL  NUMBER  OF  BRANCHES  TAKEN  =  174 

OPTIMAL  INTEGER  SOLUTION 
SUBPROBLEM  1 

18  25  6  23  16  43  0  0  0  1  1 

SUBPROBLEM  2 

39  20  20  010  101  l 

SUBPROBLEM  3 

20  1  24  23  13  20  1  0  1  0  0 

SUBPROBLEM  4 

28  12  0  2  21  4  0  1  1  0  0 

MAX  OBJECTIVE  VALUE  =  105 


i  r.  n  >*i  n  r.  u  o  n  n  a  o  c  .t  n  r.  ►-  ;i  o  o  n  o  o  o  o  n  o 
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STAIRCASE- STRUCTURED  MATRIX  RR ANCH- AND- BOUND  CODE 
PURE-INTEGER  LINEAR  PROGRAMMING  IN  GENERAL  INTEGER  VARIABLES 
SOA:  TWO  BOUNDING  PROCEDURES/  EXTRA  TIMERS  IF  NEEDED 
WRITTEN  BY  LYNNE  POLLENZ/  LAST  UPDATED  MAY  1979 

IMPLICIT  REALM  (A/C/E-H/O/P/R-W/Z)/  REALM  (B,D,X,Y), 

1  INTECERM  ( I-N/Q) 

COMMON/GESTLT/PR ICE( 130)/ ICURX(130)/ISUMC/ ITS I NY/ LI STL/ NBRANC/ 

1  US? NP/ JFCOL(ll)/ JPROW(ll )/ JFELE  M( 1 1 )/MAXC(10)/ MAXC2( 10) 

COMMON/ T I  ME RS/ 1TOT/TSTORE/TIMELP/TIMEDR/TIMEDC/TIMINV 

"AIN  PROGRAM 

START  TIMER 

00  ITOT  =  IHPTIM(l) 

INPUT  PROBLEM  DATA 

CALL  INPUT! IFPKOB/INITBD) 

IT  I MIN  =  IHPTIM(l) 

I F  ( IFPROB  •  f.Q •  0)  GO  TO  1000 
CALL  BOUNDR 
ITIMLP  =  IHPTIM(I) 

TSTORE  =  0. 

TIMELP  -  0. 

UNDER  THE  SCORE  COMPILER/  'D'  IN  THE  FIRST  COLUMN  IS  READ  AS  'C' 
UNLESS  A  SPECIAL  OPTION  IS  USED/  IN  WHICH  CASE  IT  IS  TREATED  IS  A  * 
TTMINV  =  0. 

TIMEDC  =  0. 

TIMEDR  a  0. 

NBRANC  *  0 
NtrOPWD  a  1 

APPLY  BRANCH- AND- BOUND  SEARCH  ROUTINE/  STARTING  IN  PERIOD  1 
VS  =  1 

CALL  RESTOH(l) t 
CALL  3ANDP( INIT3D) 

0  f s  (LISTL.EQ.O)  GO  TO  500 

A  PARTIAL  SOLUTION  HAS  BEEN  FOUND.  GO  ON  TO  THE  NEXT  PERIOD. 

NS  =  VS  ♦  1 
•I*OPWD  a  NFORWU  ♦  1 
CALL  RESTOK(l) 

CALL  NHFNTR 
GO  TO  30 

:,tV>  TIMER/  ALL  DONE.  REPORT  TIMING  INFORMATION. 

J  TI  v2  =  IHP*I*»{1) 

r  it  a  uri''2-tTnT)/iooooo. 

'"IT*’  (21/1)  TOT 

-1*:vaT  (/'  f(.T Al.  SJLUTIJN  TIME  =  '/Fd.2/'  SECONDS') 

W  «  r NP  =  (  ! Tlu !  V  -  T  TOT )/ 1  00000. 

-MTr.  (21/2)  Ti  Ml  NP 

*  V^Tt'  Tf  i.  r’ OP  INPUT  a'/F7.2/'  SECONDS') 

TT'iPnC  =  (ITI"LP  -  I T IMTN )/ 1000 00 . 

*  it.:  (21/3)  *imklp 

3  T 1 MF  '■')«  LP  SOLUTIONS  s'/F7.2/'  SECONUS') 

*  (Ti  (21/N)  TI*PKC 


n  o  n  r. 


20 


4  I-  1 K V A T {  *  TIME  FOR  BOUNDING  4  PRICING  s'/F7.2/'  SECONDS') 

TjTQR5,  =  TS’URF,  -  TIMFLP 

iiRlTK  (21/b)  TSTORE 

5  FORMAT ( *  Tl^t  FOR  BOOKKEEPING  OPERATIONS  =  '/F7.2/'  SECONDS') 

/SITE  (21/7)  NPORMD 

7  FT* MAT( '  NUMBER  OF  FGPMAKD  STEPS  TAKEN  =  '/l7) 

*»ITK  (21/8)  NbRANC 

FnkMAT('  TOTAL  NUMBER  OF  BRANCHES  TAKEN  ='/I10) 
n  A  R 1 T  E  (21/ V)  TIVINV 

D9  F.lRv*?('  TIME  SPENT  INVERTING  RASIS  ='/F7.2/'  SECONDS') 

C  *PITE(2t/lC)  TIMEDC 

no  FOR  MAT  ( '  TIME  SPENT  IM  SUBROUTINE  DCHUZC  ='/F7.2/'  SECONDS') 

0  4C ITE( 21/11)  TIMfcDR 

i  ll  r‘OPVAT( '  TIME  SPENT  IN  SUBROUTINE  DCHUZR  ='/F7.2/'  SECONDS') 

C  JUTPUT  OPTIMAL  SOLUTION 
CALL  WRAPUP 
GO  TO  100 
1000  STOP 
END 

BLOCK  DATA 
C 

INITIALIZE  GLOBAL  PROGRAM  CONSTANTS 

SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  BRITTEN 
BY  J.A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

IMPLICIT  REALM  (  A/C/E-H/Q/P/ R-W/7.)/  REALM  (B/D/X/Y), 

1  INTSGERM  (I-N/Q) 

COMMON/TOTSIZ/  M AXTRW/ MAXTCL/ MAXNP/M AXELE/M AXAEL 

COMMJN/CONSTS/ZTOLZE/ZTOLP  V/  ZTCOST/ 7.TOLSM/NEGI  NF/NEMAX/  NRMAX/QBL/ 
1  NLES/NTMAX/OA/QI/QF/QN/QSUB/QB/OC/QE/QH/QL/QO/QR/QM/QG 
DATA  MAXTRB/70//MAXTCL/130//M AXNP/10//MAXELE/2000//MAXAEL/1000/ 
DATA  ZTOLZE/1. E-5//ZTQLPV/ l.E-4// ZTCOST/ l*E-3/#  ZTOLSH/1. E-10/ 

DATA  NKMAX/60//NTMAX/500//NEGINFA-100000/ 

DATA  DHL/ '  '//QA/'A  '//QI/'I  '//QF/'F  '//QN/'N  '// 

1  QSUB/'SUB  '//QB/'B  V/QC/'C  V/QE/'E  '//QH/'H  '// 

?  QL/'L  '//QO/'O  '//QR/'R  '//OM/'M  V/QG/'C  '/ 

END 

SUBROUTINE  INPUT( IFPROB, INITBD) 

r 

C  READ  ALL  PROBLEM  DATA 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 

C  BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C  ** PRESCRIPT ION  OF  PARAMETERS*** 

C  TFPROfa  s  NONZERO  PROBLEM  ID  NUMBER  (OUTPUT) 

C  INITBO  s  INITIAL  LOWER  BOUND  ESTIMATE  FOR  MAXIMAL  OBJECTIVE 

C  VALUE  (OUTPUT) 

C  DATA  MUST  BE  INPUT  IN  THE  FOLLOWING  ORDER:  PROBLEM  CONSTANTS/ 

C  FOP  EACH  PERIOD:  KOWS/OPFDIAG.  COLS/DIAGONAL  COLS/RHS/BOONDS 

C  THE  OBJECTIVE  ROW  MUST  BE  THE  FIRST  ROW  IN  EACH  SUBPROBLEM 

C  VARIABLES  LOCAL  TO  ORE  PERIOD  MUST  APPEAR  AFTER  VARIABLES 

C  Of  THAT  PERIOD  WHICH  ALSO  APPEAR  IW  THE  FOLLOWING  PERIOD. 

C 

IMPLICIT  REAL*4  (A/C/F-H/O/P/R-W/Z)/  REAL*8  (B/D/X/T)/ 

1  INTEGERM  (I-N/Q) 

INTEGER  JH/K INB AS/LA/LE/ I  A/ IE 
INTEGER  NAME (6)/ NAMCOL( 130/2) 

DOUBLE  PRECISION  E(2000)/ATEMP1/ATEHP2 
MEAL  A(1000) 


nfifinorinnnncnfinomnnnnnnnnnorr.iinnnnnn 
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1 

7000 


COMMON/TOTSI 2/  MAXTRW/HAXTCL/MAXNP/M AXELE/M AX AEL 

COMMON/CONSTS/  ZTOLZE/ZTOLPV/ZTCOST/ZTOLSM/NEGINF/NEMAX/NRMAX/QBL/ 

l  nles,ntmax/QA,gi,qF/Qn,qsob,qb,qc,qe:,qh,ql,qo,qr,qm,qg 

CflMMOX  XLB(122)/XUB( 122)/DE/DP/B(60)/X(60)/ Y(60)/YTEMP(60)/A/ 

1  E/MSTAT/ IOBJ/IRONP/ITCNT/INVFRQ/ ITRFRQ/ JCOLP/ 

2  NROW/NCOL/NELEM/NETA/ NLELEM/NLET A/NULLEM/NUETA/ JH(60)/ 

3  KIN8AS(122),LA(122),LE(502),IA(1000),IE(2000) 

COMMON/GESTLT/PRICE ( 130)/ ICURX( 130)/ ISUMC/ITSINV/LISTL/  NBRAMC# 

1  NS/  MP,  JFCOLdl  )/ JFROW(ll)/ JPELEM(11)/MAXC(10)/MAXC2(10) 
COMMON/OFFD AG/OFFOC 2000 )/LCOLOO( 130)/ IR0N0D(2000)/C0STO30) 

♦♦♦DESCRIPTIONS  OF  LOCAL  VARIABLES*** 

JCPTR  POINTS  TO  CURRENT  OFFDI AGONAL  INPUT  COLUMN 

NAMCOLd/*)  =  NAME  OF  COLUMN  I 

ICS1/ICS2  =  NAME  OF  CURRENT  COLUMN 

NTCOL  =  TOTAL  NUMBER  OF  COLUMNS  SC  FAR 

NCOL1  =  FIRST  COLUMN  OF  CURRENT  SUBPROBLEM 

NRONL  =  LAST  SLACK  COLUMN  (AND  LAST  RON)  OF  CURRENT  SUBPROBLEM 
MAXTRtf  =  MAX  TOTAL  NUMBER  OF  RONS  (INCLUDING  NS  OBJECTIVE  RONS) 
MAXTCL  =  MAX  TOTAL  NUMBER  OF  COLUMNS  INCLUDING  SLACK  COLUMNS 
MAXNP  =  MAXIMUM  NUMBER  OF  SUBPROBLEMS 

MAXAEL  =  MAX.  NUMBER  OF  NONZERO  ELEMENTS  OF  ANY  DIAGONAL  BLOCK 
MAXELE  =  MAXIMUM  SUM  OF  NONZERO  ENTRIES  IN  ALL  DIAGONAL  BLOCKS 

*** DESCRIPTIONS  OF  SOME  IMPORTANT  VARIABLES  IN  BLANK  COMMON*** 

B( I)  =  RIGHT  HAND  SIDE  OF  ROW  I  C1N  CURRENT  SUBPROBLEMJ 
X(I)  =  LP  VALUE  FOR  JH(I)/  NHICH  IS  THE  VARIABLE  BASIC  IN  RON  I 
A  CONTAINS  THE  NONZERO  ELEMENTS  OF  THE  CONSTRAINT  MATRIX/  INCL. 
THE  OBJECTIVE  RON  IOBJ.  LA(J>  =  LOCATION  IN  A  OF  THE  FIRST 
ELEMENT  OF  COL  J.  IA(I)  =  RON  IN  NHICH  ELEMENT  I  OF  A  BELONGS. 
E  CONTAINS  THE  NONZERO  ELEMENTS  OF  THE  CURRENT  LP  BASIS  INVERSE 
IN  ETA  VECTOR  FORM.  LE,  IE  ARE  TO  E  AS  LA/  IA  ARE  TO  A. 

MSTAT  FLAGS  FEASIBILITY  OF  CURRENT  LP 

ITCNT  =  NO.  OF  SIMPLEX  ITERATIONS  SO  FAR;  IF  >  ITRFRQ,  STOP 
INVFRQ  =  NOMBER  OF  SIMPLEX  ITERATIONS  BEFORE  E  IS  REINfERTED 
NROW  =  NO.  OF  ROWS)  NCOL  =  NO.  OF  COLUMNS  IN  CURRENT  SUBPROSLEM 
KINSAS(J)  =  (  I  IF  J  IS  BASIC  IN  RON  1,  I.E.  J  =  JH(I) 

0  IF  J  IS  NONBASIC  AT  ITS  LOWER  BOUND  XLB(J) 

-l  IF  J  IS  NONBASIC  AT  ITS  UPPER  BOUND  XUB(J)  ) 

♦♦•DESCRIPTIONS  OF  SOME  VARIABLES  IN  GESTALT  COMMON*** 

JFCOL(I)/JFRQW(I)/JFELEM(I)  ARE  FIRST  COL/  RON/  ELEMENT  OF  k, 
RESPECTIVELY/  FOR  PERIOD  I 

NP  =  TOTAL  NUMBER  OF  PERIODS  (OR  EQUIVALENTLY/  SUBPROBLEMS) 

NS  =  CURRENT  SUBPROELSM  NUMBER 

INITIALIZATIONS 

ITCNT  =  0 
N3  =  0 
NOt’LE«  =  0 
NTCOL  =  0 
J*COL( 1 )  =  I 
JPRON(l)  =  1 
JFKLEM(l)  =  1 
DO  1  I=1/MAXTCL 
COST ( I )  =  0. 

READ  ( 20/ 7 000/ EN 0=9999)  I FPPOB/NP/ I OBJ/ INVFRQ/ ITRFRQ/ IN ITBD 
FORMAT  (4I4/I5/I10) 

IF  (IFPROU  .EG.  0)  RETURN 
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IP  (NP.GT.MAXNP)  GO  TO  9996 
NBMAX  =  4000/NP 
NLES  =  1000/NP 
IP  ( IOBJ  .EQ.  0) IOBJ  =  1 
IP  (IOBJ.NE.l)  GO  TO  9995 
IP  (INVFRQ  .EQ.  0 )INVFRQ  =  99999 
IF  (ITRFRQ  .EQ.  0)ITRFRQ  =  999999 
WRITE! 21,8010)  IFPROB 
8010  F0RMAT!/19X,  'PROBLEM  ',14) 

C 

C  INITIALIZE  FOR  READING  EACH  SOBPROBLEM 
C 

5  IF  (NS.EQ.NP)  GO  TO  25 

NROW  =  0 
ICS1  =  0 
TCS2  =  0 
DO  10  I=1,NPMAX 
10  B(  I )  =  0. 

NS  =  NS  ♦  1 
NCOL1  =  JFCOL(NS) 

READ!20,99)  (NAME( I )/ 1=1 , 3 ) 

99  FORMAT!2A4,I4) 

IP  !NAME!2 ) .NE.QSUB  .OR.  NAME(3).NE.NS  )  GO  TO  9998 

C 

C  READ  IN  DATA  FOR  SUBPROBLFM  NS 

C 

25  RRAD< 20,101)  K1,K2,K3,K4, (NAME( I ), 1=1,4) , ATEMP1,NAME!5),NAME(6), 

1  ATEMP2 

101  FORMAT! 4A1,2A4,2X,2A4,2X,F12.4,3X,2A4,2X,F12. 4) 

IF(K1  .EQ.  QBL)  GO  TO  50 
IF(K1  .EQ.  QR  .AND.  K2  .EQ.  QO)  L=1 

IF(K1  .EQ.  QR  .AND.  K2  .EQ.  QO)  GO  TO  25 

IF(K1  .EQ.  QC)  GO  TO  150 
IP(K1  .EQ.  QO)  GO  TO  160 
IP(K1  .EQ.  QR  .AND.  K2  .EQ.  QH)  L=4 

IF!  K1  .EQ.  QR  .AND.  R2  .EQ.  QH)  GO  TO  25 

1P(K1  .EQ.  QB)  GO  TO  600 
IF(K1  .EQ.  QE)  GO  TO  700 
WRITE! 21,8020)  K1,K2,K3,K4, NS 
B020  FORMAT! '  IMPROPER  INPOT  ',4A1,'  IN  SUBPROBLEM', 1 1) 

GO  TO  9999 

50  GO  TO ! 210, 320, 400, 500), L 

C 

150  L  =  2 

NR OWL  =  NCOL1  ♦  NROW  -  1 
GO  TO  25 
C 

160  L  =  3 

ICS1  =  -9999 
GO  TO  25 
C 

C  READ  ROW  NAMES  {  =  SLACK  COLUMN  NAMES) 

C 

210  NROMsNROW+1 

NCOL=NROW 
NTCOL  =  NTCOL  ♦  1 
NAMCOL!NTCOL,l)  =  NAME! 1) 

NAMCOL!NTCOL,2)  *  NAME!2) 

C 

C  TEST  ROW  TYPet  <=,»,>=,  OR  OBJ.  ROW)  SIT  SLACK  BOUNDS 
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IF(K2.EQ.QL 

.OR. 

K3.EQ.QL) 

GO 

TO 

220 

IPOC2.EQ.QE 

.OR. 

K3.EQ.Qg) 

GO 

TO 

230 

IF(K2.EQ.QG 

.OR. 

K3.EQ.QG) 

GO 

TO 

240 

IF( K2. EQ. QN 

•  OR. 

K3.EQ.QN) 

GO 

TO 

250 

GO  TO  230 

XLB(NROW)  = 

0. 

XUB(NROW)  = 

1.  E4 

GO  TO  250 

230  XLB(NROW)  =  0. 

XUB(NROW)  =  0. 

GO  TO  250 

240  XLB(NRUW)  =  0. 

XUB(NROW)  =  1.E4 
A(NROW)  =  -1. 

GO  TO  260 
250  A(NRUW)  =  1. 

26C  I  4 ( NROM )  =  NROM 

L4(NP.OW)  =  NROW 
NELEM=NRON 
GO  TO  25 
C 

C  MATRIX  ELEMENTS 

C 

320  J  =  3 
K  =  4 

IE  (DABS(ATEMPl)  .GT.  ZTOLZE)  GO  TO  324 

Js5 

K=6 

IF(DABS(ATEKP2)  .LE.  ZTOLZE)  GO  TO  25 
ATEMP 1=ATEMP2 

C  TEST  FOR  COLUMN  MATCH 

324  IF  (NAME(l)  .EQ.  ICS1  .AND.  NAM£( 2)  .EQ.  ICS2)  GO  TO  330 
NCOL  =  NCOL  ♦  1 
NTCOL  =  NTCOL  ♦  1 
ICS1  =  NAME(l) 

ICS2  =  NAME (2) 

NAMCOL( NTCOL/ 1)  s  ICS1 
NAMCOL(NTCOL/2)  =  ICS2 
LA(NCOL)  =  NELEM  ♦  l 
C  RECORD  OBJECTIVE  VALUE 

330  IF  (NAME(J)  .EQ.  NAMCOL(NCOLl-HOEJ-l,l )  )  COST(NTCOL)  =  ATEMP1 
C  TEST  FOR  RON  MATCH 

DO  340  I  =  NCOL1/NROWL 

IF  (NAME(J).NE.NAMCOL(I/l).OR.NAME(K).NE.NAMCOL(I/2))GO  TO  340 

NELEM  =  NELEM  ♦  1 

IA(NPLEM)  =  t  -  NCOL1  ♦  1 

A( NELEM)  =  ATEMPl 

LA(NCOL+l)sNELEM*l 

IF(K  .GT.  5)  GO  TO  25 

IP(DABS(ATEMP2)  .LE.  ZTOLZE)  GO  TO  25 

J  =  5 

K  =  6 

ATEMPl  *  ATEMP2 
GO  TO  330 
340  CONTINUE 

•RITE (21,8300)  NAME(J),NAME(K),NAME(1),NAME(2) 

B300  F0RMATO7HN0  MATCH  FOR  RON  ,2A4,11H  AT  COLUMN  ,2A4) 

GO  TO  9999 


C  i?Fn I agonal  matrix  elements 

i: 

400  J  -  3 
K  =  4 

IF  (DA8S( ATEMP1)  .GT.  ZTOLZE)  GO  TO  420 
J  =  5 
K  =  6 

IF(l)AnS(ATE>»P?)  .LE.  ZTOLZE)  GO  TO  25 
ATEMP1  =  4TEMP2 

C  TEST  FOR  COLUMN  MATCH 

420  IF  (NAMF(l).EQ.ICSl  • AND*  NAME(2).EQ.ICS2)  GO  TO  450 
IOSl  =  MAME(l) 

ICS2  =  NAME  (2) 

I  BEG  =  JCPTR  ♦  1 
I  END  =  JFCOL(NS)  -  1 
IF  (1PEG.GT. IEND)  GO  TO  435 
DO  430  1  =  1  BEG/ 1END 

LCOLOL  ( I )  =  NOELEM  ♦  1 

IF  ( ICS 1 . EQ. NAMCOL ( I , 1 )  .AND.  ICS2.EQ. HAMCOL( 1/2))  GO  TO  440 

430  CONTINUE 

435  !a9ITK(21,825C)  ICS1,1CS2,NS 

8250  F1kMAT('  NO  MATCH  FOR  COLUMN  *,2A4,*’  IN  SOBPROBLEM', 14) 
f.n  TO  9999 
440  JCPTR  =  I 
C  TEST  FOR  ROW  MATCH 

450  OH  460  I=NCOLl, NTCOL 

IF(NAME(J).NE.NAMCOL(I,l).OR.NAME(K).NE.NAMCOL(I,2))  GO  TO  460 

NOELEM  *  NOELEM  ♦  1 

TROUGH  NOELEM)  =  I  -  NCOL1  ♦  1 

OFFD(NOELEM)  =  ATEMP1 

LCOLOD( JCPTR  ♦  1)  =  NOELEM  ♦  I 

IF  (K.GT.b)  GO  TO  25 

IF(DAfcS(ATEMP2)  .LE.  ZTOLZE)  GO  TO  25 
J  =  5 
K  =  f 

ATEMP1  =  ATEMP2 
CO  TO  450 
ioO  CONTINUE 

«^ITE( 21/8300)  NAME(J),NAME(K),NAME(1),NAME(2) 

Gn  TO  9999 
C 

l  :<IGilT  HAND  S1DF. 

C 

600  J  =  3 
K  =  4 

IF  (DABS(ATFMPl)  .GT.  ZTOLZE)  GO  TO  530 

J=5 

K=6 

IF( DABS( ATEMP2)  .LE.  ZTOLZE)  GO  TO  2b 
ATEMP1=ATEMP2 

C  TEST  FOR  ROW  MATCH 

630  DO  540  I  =  NCOL1 /NROWL 

IF<NAMF(J).NE.NAMCOL( 1,1 ) . OP. NAME(K) . NE. NAMCOL(I,2))GO  TO  540 

B(I-NC0L1*1)  =  ATEMP1 

IF(K  .GT.  5)  GO  TO  25 

IF(DABS(ATErfP2)  .LE.  ZTOLZE)  GO  TO  25 

J  =  5 

K  =  6 

ATEMP1  =  ATEMP2 
GO  TO  530 
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MO  CONTINUE 

ViRITE(  21,8300)  NAM£( J),NAME(K) 

GO  T!)  9900 
C 

C  BOUNDS  ON  INTEGER  VARIABLES 

C 

*00  K  =  NRUW  ♦  1 

C  TNPUT  LOWER  AND  UPPER  BOUNDS  ON  DECISION  VARIABLES 

READ  (20,3600)  (XLB(J)/  J=K/NCOL) 

RTAI)  (20,8500)  (XUR(J),  J=K,NCUL) 

*500  FORMAT  (15FE.0) 

C  CHECK  PROBLEM  SIZE 

JFCULUS+l)  =  JFCOL(NS)  ♦  NCOL 
JFK0W(.NS4l)  =  JFROW(NS)  ♦  NRON 
JEELEMOlMl)  =  JFELEM(NS)  ♦  NELF.M 
IF  (NROV.Gl.MRMAX)  GC  TO  9996 
Ir  (9KLKM.GT.“AXAEL)  GU  TO  9906 
IF  <JFUn,(NMl)  .GT.  (MAXTRW  ♦  1))  GO  TD  9996 

!=■  <JFCnc(NS+l)  .GT.  (MAXTCL  ♦  1))  GO  TO  9996 

IF  (JFELt-.'*(vs*l)  .GT.  (MAXELL  «•  1))  GO  TO  9996 
JC  1JT!<  =  JFCHL(NS)  ♦  NROW  -  1 
C  *  LL  I‘iPGTO 
GO  TO  5 

\ 

r  •  n  nF  INPUT 
C 

700  in. ? r1  U  i  =  JFPnW(NP«-l)  «  1  -  *P 

VJCOLS  =  JFCOL( NP+1)  -  JFROW( NP+1 ) 

'*..'<IT  :(  21/8O00)  NTRON/NSCOLS/NP 

RHOP  rr],<MAT(5X/l<!/'  ROWS 2X, 1 4, '  COLUMNS  *,2X, 14, *  PERIODS*//) 

txFTMRN 

C 

C  ERROR  MESSAGES  FOR  INPUT  ERRORS 

5,' 

999r  *9 1  TEC  21/0  960) 

°0t> C  p'1RMAT('  UJECT1V-  ROW  MUST  HE  THE  FIRST  ROW*) 

GO  TO  9999 

9996  \y I TE ( 21/ 0970) 

F97't  r  UHATC  *  PROBLEM  IS  TOO  LARGE  FOR  CURRENT  DIMENSIONING*) 

*1-*  IT*. <21/ 8975)  MAXNP/MAXTKW/M  AXTCL/M  AXELE 
"976  FOP MT(  ’MAX  SU8S  =  *,  13/  'MAX  ROWS  =  '/ 14/ 'M AX  COLS=  */ 14/ 'MAX  ELE=',I6) 

r,  »  XO  *5^99 

>99  .  rtCITKC  21/3990)  NS 

y.  9  90  F1PWATC  SUi'rPUBLFM*/  14/  *  MJT  IN  PROPER  POSITION*) 

9 r>  *  IPRRfJP  -  0 
RETURN 

S’*'JPOUTI\E  INPSTO 

C  STORE  ALL  DATA  KrLFVANT  TO  SUBPROULEM  NS 

C 

IMPLICIT  REALM  (A/C/F-H/O/P/R-W/Z)/  Rt AL*8  (B/D/X/Y)/ 

1  INTEGER**  (I-'I/D) 

I  NTEGER  JII/V  1NBAG/LA/LF./IA/IE /LCDLA, IRUW A, IESTOR/LESTOR 
DOUBLE  PRECISION  S< 2000  ) / ESTOP e, ( 4 C 00 ) 
r "AL  1(1000) 

CO" '  JN/CONSTS/  ZTOLZE/ZTOLPV/ZTCUST/ZTULSM/NECINF/NEMAX/NRMAX/QBL/ 
I  NLr'S/NT.NA  X/  QA, «I /  OF,  ON,  QS Ul>/  Wli,  >iC,UF ,GK,UL,GO,QR,GM,WG 
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COMMON  XLRd22),XUB(122),DE,DP,n(60),X(60),Y(60),mNP(50),4, 

1  t,MS  TAT/  IOPJ, TkOMP,  ITCNT, INVFRQ,  ITRFRQ, JCOLP, 

2  N ROW, NCOL, SKLEM,NET A,NLEL£M,NLETA,NUELEM,NUETA, JH(60), 

3  XINPAS< 122 ),LA(122),L5( 502), 1A(1000),IE( 2000) 

CQkMON/G:  SILT/ PR ICE ( 1 30 ),  ICURX( 1 30 ) , 1SUMC, ITS  I W,LISTL, NBR ANC, 
1  NS,NP,JFCnLdl),JFkOhdl),JFELEM(ll),MAXC(10),MAXC2(10) 

C  1MMCJN/T  LMZRb/lT0T,T5T0RK,TT  HELP,  TIHEDR,TIMEDC,T  ININ  V 
COkMUN/SDBSAV/BMUD{70  ),  90RIG(7O ),  X STORE  (  70),XU6STO(1 30  ), 

1  X!  BSTCK  130),  X'JBORG  ( 130  ),XLBOPG  (  130),  AS  TORE  (  2000),  E  STORE, 

2  T POVA( 2000), LCOLA(l 40 ),IKSTOH( 4000), LESTORC 1002), JHSTOR(70), 

3  I NB ST 0(130),  IP  ARTC(  10  ),INVSTG(10  ),  NELSTOdO  ),  NETSTO(  10) 


C  C n M M 1 N  BLOCK  SUrSAV  IS  A  storage  arf.a  for  the  subpkqblem  variables 

r.  XL U'JKG,  XU  PllRG, PORI  G  STOKE  THE  ORIGINAL  VALUES  OF  THE  CORRESP.  VARS. 

C  .  -  rj  0(1)  =  RHS  FUR  ttO 1  I  AFTER  QFFP I  AGONAL  VALUES  HAVE  BEEN  ADDED 


c 

C  AFTER  READING  INPUT,  STOKE  THE  DIAGONAL  MATRIX  ELEMENTS  AND  RHS 

C 

IPTR  =  JFELEM(NS) 

LENGTH  =  JFFLSH(NS  ♦  1 )  -  IPTR 
DO  10  1=1, LENGTH 

ASTORL(IPTR)  =  A(I) 

IKONA(IPTR)  =  I A( I ) 

10  IPTR  *  IPTR  ♦  1 

1 PTK  =  JFCOL(NS)  -  2  ♦  NS 
I END  =  NCOL  ♦  1 
DO  15  1=1, IEND 

15  LCOL A( IPTR  ♦  I)  =  LA(I) 

C  STORE  ORIGINAL  RHS  B;  BOUNDS  XUB,XLB 
I DT  R  =  JFROM(NS) 

LO  20  1  =  1,  NROVf 

30RIC ( IPTR)  =  Bd) 

20  IPTR  =  IPTR  ♦  1 

I°TP  =  JFCOL(NS) 

DO  30  1=1, NCOL 

XUBORG(IPTR)  =  XUBd) 

XLBORC(IPTR)  =  XLB(I) 

30  IPTR  =  IPTR  ♦  1 

RETURN 
C 

C  STORK  THE  CURRENT  STATE  OF  SUBPROBLEM  NS 

C 

KNTRY  STORE 

C  STORE  THE  PARTIAL  OBJ.  VALUE  ISUHC;  RHS  b)  SOLUTION  X 
IT1ME  =  IHPTIM(l) 

IPAPTC(NS)  =  ISUMC 
I°TR  =  JFROW(NS) 

DO  300  1=1, MROW 

BMOD(IPTR)  =  6(1) 

XSTORE(IPTR)  =  X(I) 

JHSTOR(IPTR)  =  JH(I) 

300  IPTR  =  I PTK  ♦  1 

C  STORE  BASIC  VAR  INDICATOR  KINHASJ  BOUNDS  XUR,XLB 
IPTR  =  JFCOL(NS) 

DO  350  1=1, NCOL 

INBSTO(IPTR)  =  KINBAS(T) 

XUPSTOdPTR)  =  XUQ(l) 

XLBSTO(IPTR)  =  XLB(I) 

350  IPTR  =  I FTP  ♦  1 

C  STORE  BASIS  INVERSE.  IF  TOO  LANGE,  REINVENT. 
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IF  (NETA.LT. NLES )  GJ  TO  390 
CALL  INVERT 
ITSTNV  =  0 

390  INVSTO(NS)  =  ITS INV 
NELSTO(NS)  =  NELEM 
NETSTO(NS )  =  NETA 
LSPfk  =  (NS-1)  *  NLES 
JEPTR  -  (NS-1)  *  NSMAX 
00  400  1=1, NELEM 

JEPTR  =  JEPTR  ♦  1 
KSTURE(JEPTR)  =  E(I) 

400  lESrOR(JFPTR)  =  IE( I) 

ISND  =  NETA  ♦  1 
DO  450  1  =  1/  TEND 

LKPTR  =  LEPTR  ♦  1 
450  LESTOP(LEPTK)  =  LE(I) 

IdMF?  =  IHPTIMd) 

TSTOFE  =  TSTORh  ♦  ( ITIME2-ITIME)/I0000G. 

RETURN 

^ND 

SUBROUTINE  RESTOR(MNFLAG) 

C  RESTORE  ALL  DATA  KELEVAMT  TO  SUBPROBLSM  NS 

C  ***DFSCKIPTION  0»  PARAMETER*** 

C  MNFLAG  =  0  IF  RACKTRACKINC,  1  IF  TAKING  A  FORWARD  STEP  (INPUT) 

r  IF  TAKING  A  FORWARD  STPP,  RESTORE  ORIGINAL  LP-OPT.  IF  TAKING 

C  THIS  ST^P  FOR  THE  FIRST  TIME/  COMPUTE  ORIGINAL  LP  SOLUTION. 

r 

IMPLICIT  REAL* 4  ( A/C/ E-H,a,P,R-W, Z),  REAL*8  <8,0, X,Y), 

1  lNTFGe‘k*4  (I-N/Q) 

INTEGER  JH/KI’JBAS/LA/LE/  I  A/  IS/LCOLA/  IRON  A/ ISSTOR/LESTOR 
P  JUPLE  PRECISION  E(  2000 ),ESTORF.( 4009 ) 

R  •’ A L  A(1000) 
r 

COMPUN/CCNSTS/  7.T0LZE,ZT0LPV,7TC0ST,ZT0LSM,  NEGINF,NENAX,NRNAX,QBL, 
1  NLES/NTMAX/QA/Ol/OF/QN/ QSUB/ QE/ QC/QE/QH/QL/ QO/QR/QM/OG 
COMMON  XLR( 1 22), XU8< 122), DE, DP, R{ 60), X(60),Y(60),YTEMP(60), A, 

1  E/MSTAT/  KJBJ, TKOWP,ITCNT,INVFRJ, ITNFRQ, JCOLP, 

2  NROW,NCOL, NELEM, NETA, NLRLEM,NLETA,NUELEM,NUETA,JH(60), 

1  K  INBASd  22  )/LA(122)/LE(502)/IA(  1000  )/IE(  2000) 

CI*M'i)M/GKSTLT/PRICE(130),IC'JKXd30),TSUMC,ITSINT,LISTL,NBRANC, 

1  US, Np/JFCOL(l 1 )/JFKOW(ll )/JFELEM(lI)/MAXC(10 )/MAXC2(10) 

C  IMMON/TIMF.I  S/ITOT/TSTORE/TIMELP/TIMEDK/TIMEDC/TIMINV 
C  l!:.^nN/SUBSAV/bMUD(70)/bORIG(70)/XSTORE(70)/XUBSTO(130)/ 

1  XLf)STO(  1  30  )/XUPOPG(  139 ),XLBORO  ( 1 30  ),  AS  TORE  (2000  )/ESTJRE, 

2  IwOM^(20CO)/LCOLA(140)/IESTOR( 4000 ),LESTOR( 1002 ),JHSTOR (70), 
o  INBSTOd 30)/IPARTC(10)/TNYSTO(10)/NEL3TO(10)/NETSTO(10) 

C 

L  'ESTOPt  SUrtPrOHLtM  DIMENSIONS)  A  MATRIX 
r 

ITIMF  =  1  HPT  IP ( 1  ) 

r.PUW  =  JFPQW(NSKl)  -  JFROF(YS) 

KCOL  =  JFCOL(MSd)  -  JFCOL(NS) 

1 PTR  =  JFELrP( NS) 

LENGTH  =  JFFLEM(NSd)  -  IPTR 
C  RESTORE  F.LffMEMS  Of  A  MATRIX  FROM  DIAGONAL  BLOCK  NS 
i)fJ  10  1  =  1,LFKGTH 

A(T)  =  ASTORE(IPTK) 

I  A(  I  )  =  TRQhAdfTf) 
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10  IPTR  =  I P1H  ♦  1 

ipth  =  jrcoi-(vs)  ♦  ns  -  * 
i?;;n  =  r<ic dl  ♦  l 

D)  15  I=1,IFND 

15  LA(  I  )  =  LCQLAI  IPTR  ♦  I) 

JFPTk  =  (NS-1)  *  NF.MAX 
LEPTk  =  (YS-l)  *  NL3S 
MSTAT  =  URL 

IF  (HNFLAG.FQ.l)  GG  TO  200 
C 

C  BACKTRACK  STEP:  RESTORE  RHS  B;  LP  SOLOTTGN  X,JH,KINBAS;  BOUNDS 
C 

IPTR  =  JFROn(NS) 
m  20  1=1, NPO# 

'(I)  =  BMUD(tPTR) 

X ( I )  =  XSTOKK (IPTP ) 

j'Hi)  =  jhstoriiptp) 

10  IPTR  =  IPTR  ♦  1 

IPTR  =  JFCGL(NS) 

DO  40  I=l,NCOL 

xpb(I)  =  xuosrnt iptr) 

XLB(I)  =  XL)STO(IPTK) 

K l NBAS! I)  =  TNBSTQI  IPTR) 

40  IPT!'  =  IPTR  ♦  1 

c  kfstcjpe  partial  ouj.  value  if  backtracking  (not  for  forward  steps) 

ISUMC  =  IPARTC(NS) 

c 

C  RESTORE  LP  BASTS  INVERSE 

r* 

v. 

TTSiriV  =  INVSTU(NS) 

K’ELSM  =  JELSTn(NS) 

META  =  NETSTO(NS) 

DO  110  1 =l,NELEM 

JEPTB  =  JEPTR  ♦  1 
F(I)  =  ESTOKEI JEPTR) 

110  IB(I)  =  IESTOR(JEPTR) 

IFND  =  NETA  ♦  1 
DO  120  1=1, TEND 

LEPTR  =  LEPTR  ♦  1 
12C  LEI  I)  =  LESTOR(LEPTR) 

CO  TU  500 
C 

C  FORWARD  STEP 
C 

C  PFSTOHE  ORIGINAL  BOUNDS  FOR  A  FORWARD  STEP 
200  I*»TR  =  JFCOL(NS) 

00  210  I=l,*?COL 

XUB (  I  )  =  XUHORGI IPTR) 

XLB(T)  =  XLBOKGI IPTR) 

21C  IP^R  =  IPTR  ♦  1 

IPTR  =  JFROW(NS)  -  1 
DO  220  1=1, NROW 

220  D(T)  =  BOhlGI IPTR  ♦  I) 

C 

c  solve,  lp  relaxation  of  subproblem  ns 

c 

IF  (NS.GT.l)  CALL  FIXPHS 
C  LP  BASIS  STARTS  IFF  AS  ALL  SLACK  BASIS 
1TIMLP  =  IHPTIH(l) 

DO  310  I*l,NFOW  | 

. I 
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3ic  jhu)  =  r 

on  320  I=l,NkO- 
320  KINBAS(I)  =  I 

NROWP1  =  NROr  ♦  1 
DO  330  1=NRQNP1,NC0L 
33C  KIN3AS(1)  =  0 

£ 

C  SOLVE  LP 

C 

ITSINV  =  99999 
C ALL  NORMAL(ITSINV) 

IF  (MSTAT.EQ.QN)  GO  TO  2000 
iriML2  =  IHPTIM(l) 

Tf RSLP  =  TIMLLP  ♦  ( ITIML2-ITIMLP) /100000. 

C 

r.0C  ITIME2  =  IHPTIk(I) 

TSTOKE  =  TSTGRE  ♦  ( ITIME2-ITIME)/100000. 

kETUPN 

C 

r  L°  IS  INF? AS1°LE 

c 

2000  IF  (NS.FQ.l)  W'<ITE(21,2010) 

2010  FORMATC  SUBFROBL^  1  IS  INFEASIBLE') 

G I  TO  iOC 

End 

C - - — - - - - - — - — - - - 

subroutine  fixrrs 

c 

C  GIVEN  SETTING  OF  VARIABLES  FOR  PERIOD  (NS-1),  COMPUTE  NEW 

C  RHS  F'JK  PFKTOD  NS  IN  PREPARATION  FOR  A  FORWARD  STEP. 

C 

IMPLICIT  REAL* 4  <*,C,E-H/0,P,R-W,Z),  RLAL*8  (P,D,X,Y), 

1  1 NTEGEhM  (I-N,Q) 

INTEGFK  J  I/NINBAS/LA/LE/IA/IE 
DOUBLE  PRECISION  E(2000) 

REAL  A(1000) 

r.  vmmun  XLR(122),XUri(122)/0E,DP,R(60),X(6  0),Y(60),YTEMP(60),A, 

1  K/TSTAT/ IOBJ/ IROdPf ITCNT/ I NVFRQ# ITRFRQ/ JCOLP/ 

?  N»nW,NCOL,NELEfc,NETA,NLFLEM,NLETA/NUELEM,NUETA,JH(60)« 

3  KINBAS(122)/LA(122)/LF.(502)/IA(1090)/IE(2000) 

COMMOfi/GLSTLT/PHlCE(130)/TCURX(l30)/ISUMC#ITSTMV*LI  STL/  NBRANC/ 

2  IS,NP,JEcOL(il),JPROmi),  JFELKMUl  ),  MAXC(10),MAXC2(10) 

COM**  J*!/n?FDAG/OFFD{ 2000 ),LCOLOD(  130)/ IROWOD(2009),COST(  130) 

>:  START  AT  FIRST  NONSLACK  COLUMN  OF  LAST  PfcPIOD 

J’EG  *  JFCOL(NS-l)  ♦  (JFRON(NS)  -  JFRO*( NS- 1 ) ) 

J-iO  =  JECUL(NJ)  -  1 
C 

C  ORPUTE  CONTRIBUTION  TO  RES  FROM  LAST  PERIOD 

(. 

PI  100  J=J8PG,JENr> 

C  UNPACK  nFKDT/'.L’  AL  COLUMN  J  IF  NONZERO 
IKEG  =  LCOLOD(J) 

I  END  =  LCULOP(J*l)  -  1 
IF  (  IBEvj.  GT .  IfcND  )  GJ  TO  100 
DO  20  !=l/NROW 
YTFMP(I)  =  0. 

0* .  3'  I  =  IBF.Grin;JO 
tk  =  fhOWOD(I) 


;o 


onnr.omnomn  n  n  a  n  n  r. 


jO  YTEMP(lR)  =  OFFD(  I) 

C  -i'JLMPLY  COLUMN  J  BY  ICURX(J),  f  HF  VALUE  OF  X(J)  AS  LAST  SET 
VTEMP(IOBJ)  =  0. 

xvalue  =  icunx(j) 

no  40  1=1 /NR ON 

40  H( I )  =  R(I)  -  (YTFMP(I)  *  XVALUE) 

100  CCN  flNUE 

RETURN 
END 

C - - - 

SU! ROUTINE  OOUMOR 

f  CALCULATE  2  LP  BASKO  Dm  NOS  ON  FAX  OBJECTIVE  VALUES  FOR  PERIODS 

C  NP/NI-l/.../K  FOR  K  =  HP,.. .,7  TNCLUDING  OFFDI AGONAL  COLUMNS* 

C  ALSO  PRICE  OUT  THE  OFFOI AGONAL  COLUMNS. 

r  THIS  SUBROUTINE  ASSUMES  OBJKCTIVp  RON  =  1  IN  EACH  SURPR08LEM/ 

C  AND  THAT  VARIABLES  LOCAL  TO  PFR TOD  NS  ARE  NUMBERED  AFTER 

THOSr,  VARS.  WHICH  nAVE  NONZERO  ENTRIES  IN  PERIODS  NS  A  NS+1* 
ALSO  TOTAL  NUMBER  OF  NONZERO  ELEMENTS  IN  THE  CONSTRAINT 
MATRIX  INCLUDING  OFFDT AGONAL  COLUMNS  MUST  NOT  EXCEED  1000 
(FOR  CURPENTLV  DIMENSIONING). 

IMPLICIT  REALM  (A/C/E-H/O/P/R-W/Z)/  REALM  (B,D,K,  Y), 

1  INTFGER*4  ( T -N/O ) 

INTEGER  JH/KINBAS/LA/LE/IA/ IE/LCOLA, IROWA/ I ESTOR/LESTOR 
DOUBLE  PRECISION  K( 2000 )/ FSTQPE( 4000) 

PEAL  A(1000) 

C'lMMON/caNSTS/ZTOLZE,ZTOLPV,ZTCOST,ZTOLSM,NEGINF,NENAX,  NRMAX/QBL/ 
l  NLES/NTMAX,QA,OI,9F,QN,nSUS,;iB,QC,QE,QH,QL,QO,‘lR,QM,QC 
COMMON  XLB(122>/XUB(122)/DE,OP,«(60)/X(6G)/Y(60)/YTENP(60)/A/ 

1  t'/  MSTAT/  IObJ/  I  KOWP/  ITCNT/  I  NVFRQ,  I TKFRQ/  JCOLP, 

2  NRQV,NCUL,NELEm,NETA,NLELEH,NLETA,NUELEM,NUETA,JH(60), 

1  KINBAS(122)/LA(122)/LF<‘5  02)/IA(10(>0)/IE(2000) 

C'lMMlN/Gt3TLT/PRICE(130)/ICUKX(130)/ISUHC/  ITSINV/LISTL/  NBRANC, 
l  NS/NP/JFCOL(U )/JFROW(Il )/JFELEM(tl )/MAXC(10 )/MAXC2( 1 0) 
COMMON/SinSAV/BMOD(70),HURIG(70),XSTCRE(70)/XUBSTO(130)/ 

1  KLRSTO(130)/XUBORG(130),XLBOMG(130), ASTORE( 2000)/ESTORE/ 

2  IPOWA(2000)/LCOLA(1 40 ),IESTOR( 4000 ),LESTOR( 1002 ),JHSTOR(70), 

3  INBSTO(130)/IPARTC(10)/INVSTQ( 1 0 ) /NFLSTO(10 )/NETSTO( 10 ) 

COMMON/ JFFDAG/OFFO( 2000 )/LCOLOD( 130 )/IROWOD(2000)/COST( 130) 

MXC(fC)  s  WEIGHTED  LP  BDUMU  ON  MAX  OBJ  VALUE  FOR  SUBS  K+1/.../NP 
OFFDI ACONAL  VARIABLES  OF  SU3PROBLFM  Ml  HAVE  A  WEIGHT  OF  1, 

ALL  OTHER  VARIABLES  HAV&  A  WEIGHT  OF  2.  THUS  AN  LP  BOUND  ON 
PH'-:  ObJ.  VALUE  FOR  DIAGONAL  BLOCK  K  ♦  MAXC(K)  IS  A  VALID 
30'IND  ON  DOUBLE  THE  OBJ.  VALUE  FOR  PERIODS  K/K+l/. ../NP. 

(IN  TECH.  REPORT  BY  POLLENZ  C19B0J,  THIS  CORRESPONDS  TO 
MA  (CO/2/K).  ) 

MXC2(K)  =  LP  BOUND  ON  MAX  OBJ  VALUE  FOR  SUBS  K*1,...,NP.  PERIOD  K 
VARIABLES  HAVE  COST  0,  THOSE  FROM  PERIODS  >  K  HAVE  ORIGINAL  COSTS. 

(  THIS  CORRESPONDS  TO  MAXC(0/K)  IN  THE  TECH.  REPORT.) 

MAXC(NP)  =  0 
MAXC2(NP)  =  0 
LASTC  =  JFCOL(NPwl)  -  1 
DO  20  J  =  l, LASTC 
20  PRICF(J)  =  0. 

NROW  =  0 

C  n  FA  J  IN  SLACK  COLUMNS)  SKIP  OBJ.  SLACKS  FOR  PERIODS  2/.../NP 
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Dn  40  NS=1/NP 

TBEG  =  JFELEM(NS) 

I  END  =  IBEG  ♦  JFROW(NS*l )  -  JFROW(NS)  -  1 
IF  (NS.NF.l)  IBEG  =  IBEG  ♦  1 
DU  30  I=I8EG/IEND 
NPOH  =  NRQN  ♦  1 
A(NROW)  =  ASTORS(I) 

LA(NROW)  =  NRUW 
IA(NROW)  =  HR 0* 

JH(NROM)  =  NROV 
KINBAS(NROW)  =  NKOrf 
3(NRO«)  =  0. 

JCCL  =  JFCQL(NS)  ♦  1  ♦  I  -  IBEG 
IF  (WS.SU.l)  JCOL  =  JCOL  -  l 
XUB(NRDW)  =  XUBf)RG(  JCOL) 

30  XL3( NPGW )  =  XLBORG(JCOL) 

*.0  CONTINUE 

NS  =  NP 
NCJL  =  NRIJW 
NAPTR  =  NROW  ♦  1 
C 

C  <\F  AD  RIG.IT  HAND  SIDE  FOR  PERIOD  NS 

C 

100  1FXTRA  =  NS  -  1 

1'ic.G  =  JFRGN(NS)  -  MS  ♦  2 
IFND  ^  JFROW(NS+l)  -  NS 
DO  150  IRUb=IFLG/IEND 
150  J(  T  ROW )  =  DOKIGdROW  ♦  IEXTKA) 

C 

C  ADD  GN  F'JLL  COLUMN  FOR  NONSLACK  PFRIOD  NS  VARIABLES 

C 

JREG  =  IFCOL(NS)  ♦  ( JFROW(NS*l)  -  JFROW(NS)) 

JMD  =  JFCOL (NS*  1 )  -  1 
SO  300  J=JBFG/ JFND 

r  FIRST  ADD  ON  OFFDI AGON AL  ELEMENTS  OF  COLUMN  J  IF  ANY 

NCOL  =  MCQL  ♦  1 
LA(NCOL)  =  NAPTR 
IF  (NS.EQ.NP)  GO  TO  220 
TBF.C  =  LCOLOD(J) 

TEND  =  LCQLOD( J+l)  -  1 
IF  ( TbFG.GT.IEND)  GO  TO  22C 
UJ  210  1= IBEG/ IEND 

IF  ( IROWOD(I).EQ.IOBJ)  GO  TO  210 
A(NAPTft)  =  OFFD(I) 

IA(NAPTR)  =  IROWOD(I)  ♦  JFRObUS*l)  -  (NS*1) 
NAPTR  =  NAPTR  ♦  1 
*10  CONTINUE 

C  NEXT  ADO  ON  DIAGONAL  ELEMENTS  OF  COLUMN  J 

l '’EG  =  LrOLA(J*NS-l)  ♦  J»fLLSM(NS)  -  1 
T ?•’ D  =  LOJLA(J*NS)  ♦  JFcLEV(NS)  -  2 
D  l  250  I  =  IBEG/ TEND 

A(NAPTR)  =  ASTORE(I) 

IF  (lROWA(I).NE.IOilJ)  GO  TO  230 
A(NAPTk)  =  A(NAPTP)  *  2.0 
I A (NAPTR)  =  IORJ 
GO  TU  240 

230  I A(NAPTP)  =  IRONA(I)  ♦  JFRON(NS)  -  NS 

240  NAPTR  =  NAPTR  ♦  1 

250  Cl.JTP’UE 

¥  fNB 'S(NCOL)  =  0 
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XUB(KCOL)  =  XUBGRG(J) 

XL3(NC?L)  =  XLPOPG(J) 

300  CONTINUE 

LA(  NCOLM  )  =  NAPTR 

C  RESOLVE  LP  *1  Trt  PERIOD  NS  VARIABLES  ADDED 
ITS INV  =  99999 
CALL  NURMAL(ITSINV) 

<• 

C  add  OFf DI AGON AL  COLUMNS  FROM  PERIOD  NS-1  TO  LP 

C 

NAPTHO  =  NAPTR 
NTEMPC  =  NCOL 

C  START  AT  FIRST  NQNSLACK  COLUMN  OF  PERIOD  NS-1 

J0;:c  =  JFC3LCNS-1)  ♦  (JFKQfc(NS)  -  JFROKNS-1)) 

JEUD  =  JFCOL(MS)  -  1 
DO  430  J= J3EG, JENO 

I  i!EC  =  LCOLOD(J) 

IEND  =  LCOLOD(JtI)  -  1 
NCQL  =  NCQL  ♦  1 
LA( NCOL)  =  NAPTRO 
KINBAS(NCOL)  =  0 

C  ADD  IN  COST  OF  0  AND  BOUNDS  OF  COLUMN  J 
IA(NAPTRO)  =  1 
l( NAPTRO)  =  0. 

NAPTRO  =  NAPTRO  ♦  1 
XLH(NCOL)  =  XLBORG(J) 

XUB(NCOL)  =  XUEORG(J) 

IF  (IEEG.GT.IEND)  G3  TO  430 
OU  470  1= IBEG, IEND 

IF  UPOWOD(I).EQ.IOBJ)  GO  TO  420 
IA(NAPTRO)=IRONOD(I)+JFRCM(NS)-NS 
A( NAPTRO)  =  OFFD(I) 

NAPTRO  =  NAPTRO  ♦  1 
42C  CONTINUE 

410  CONTINUE 

LA( NCOL+1 )  =  NAPTRO 
C 

C  OUT  EACH  OFFOIAGONAL  COLUMN  OF  SUB.  NS  AND  STORE  IN  PRICE. 

C 

CALL  UPDATX 
CALL  FOR. VC 
CALL  qTRAN 
DO  500  J=JBEG, JENO 
DFPICE  =  0. 

IBEG  =  LCOLOD(J) 

IEND  =  uCOLOD( J*1 }  -  1 
IF  (TBEG  .GT.  IEND )  GQ  TO  470 
DO  450  Is IBEG, IEND 
IP  s  IKOWOD(I) 

IF  (IR  .EQ.  IOBJ )  GO  TO  450 
IH  =  IK  ♦  JFkOM(NS)  -  NS 
DF  s  OFFD(I) 

DPRICE  =  DPR ICE  ♦  (DE  *  Y(IR>> 

450  CONTINUE 

C  *rC0KD  REDUCED  COSTS,  RECALLING  THAT  OBJECTIVE  VALUES  HERE  DOUBLED. 
470  PRICE(J)  =  DPRICE  /  2.0 

bOO  CONTINUE 

*RITE(  21,  lCIO  NS 

HR 1TE( 21,1020  <  PRICE!  J),J=JI>FG,JEND) 

1010  FORMAT! *  PRICK  OF  OFFDIAGONAL  COLUMNS  IN  SUB. ',13) 
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1020  P0kMAT( 10F7. 2) 

f* 

l, 

C  RESOLVE  LP  WITH  OPFDI AGONAL  COLUMNS  ADDED  ON.  CALCULATE  MAXC2. 

C 

1TSINV  =  9S999 
CALL  NIRMAL(ITSINV) 

UOIt.l  =  X(ICBJ)  ♦  ZTOL7E 
IE  (DOBJ.GE.O.)  GO  TO  510 
DnB J  =  X(TOHJ)  -  7.TCLZE 
31 0  MAXC2( NS-1 )  =  ID1NT(DOLJ/2.0) 

C 

C  AOT  IN  OFFDI  AGnJ.  AL  COSTS  FOR  CALCULATION  OF  MAXC 

C 

NCOL  =  NTEMPC 
OH  520  J=JL EG, JEND 

NCOL  =  NCOL  ♦  1 

52C  A(LACNCOL))  =  COST(J) 

C 

C  RESOLVE  LP  WITH  OFFDI AGONAL  COLUMNS  ADDED  ON  AND  COSTS  SET. 

C 

ITSINV  =  99909 
C*LL  R07MAL(ITS1NV) 

IF  (MS  rAT.E9.9N)  GJ  TO  2000 

US  =  NS  -  1 

MCQL  =  NTSMPC 

L'3'jJ  =  X(IOBJ)  ♦  ZTULZE 

MXC(NS)  =  IDlNT(DOnj) 

IF  (On3J.  LT.  0.)  .lAXC(NS)  =  MAXC(NS)  -  1 
I>  (NS.Eh.l)  GO  TO  100 

MAXC  AND  PRICE  HAVE  HEEN  CALCULATED.  RETURN. 
r 

K  =  LP  -  1 
WI?ITE(  21/  1505) 

150  F"*RMAT</#  MAXIMUM  CUMULATIVE  OBJECTIVE  VALUES') 

k.  R  ITE(  21/ 151 C  )  (MAXC(I)/I  =  1/  K) 

1510  I’ORvATC  FOP  LAMBDA  =  .5:', 916) 
flMIT7(21,15?Ci)  (MAXC2(I)/I  =  1/K) 

1520  F’JHMATC  FOP  LAMBDA  =  0:'/9I6) 

i»^ITE(21/152b) 

1525  FHPP\T('  ') 
tv  ’T'jr’V 

K 

C  SUPPMOILE.'  NS  IS  INFEASIBLE.  QUIT 

f* 

2000  W'JITE(21/2C10)  NS 

2010  FORMATS  SUPPRUHLEM'/ 13^ '  IS  INFEASIBLE') 

STOP 

FNO 

C - - 

SUBROUTINE  updatx 

C  UPDATr.  PFo  LSI  MG  VARIABLE  HOUNDS  ACCORDING  TO  YINBAS*  THEN 

L  USE  «A3IS  INVERSE  T(J  TRANSFORM.  INTO  CORRECT  X  VECTOR. 

C 

IMPLICIT  veil* 4  (A,C,F-H/0,P/R-*/7),  RE AL*3  <B,D/X/V>/. 

1  INTEGER *4  (l-v,Q> 

INTI  Gr.  -M#  F  INl’AS/LA/Lfr'/IA/IV 
D'vn L-.  P^risiON  7(2000 

h-;al  A (  1  000 ) 


L 


n  o  f;  o  ->  r. 
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C 

COM*  IN  XL’U122)/X,'H(122)/DF.,f>n/*(60)/r(6  0)/Y(60)/VTEMP(6Q)/A/ 
l  E/^STAT/TfU  J/  tkO.VP/  ITCNT/  I  KVFRQ/  ITR*KG/ JCDLP/ 

1  NRUR/NCOL/NFLEM/NETA/NLFLEF./NLrTA/NUELEM/NUETA/JH^O)/ 

1  KINRAS(122)/LA(122)/LF(502)/IA(1000)/TE( 2000) 

C 

CALL  SiIFTR(l/3) 

DO  100  J=  1/ •'ICOL 

IF  (*,INR4S(J))20/40,100 


20 

DE 

=  XUC(J) 

GO  Tp  60 

to 

D.» 

=  XLB(J) 

fjO 

LL 

=  L  A  (  J  ) 

KK  =  L A( J*1 )  -  1 
00  80  I=LL/KK 
I R  =  I  A(  T  ) 

•10  Y(TR)  =  Y(IR)  -  A( I )*Di 

10C  CONTINUE 

CALL  FTRA'Kl ) 

CALL  SHIFTKC3/2) 

RETURN 

END 

SUBROUTINE  FTRAN(IPAR) 

C 

PERFORM  FORWARD  TRANSFORMATION  ON  COLUMN  STORED  IN  VECTOR  Y 
SUBROUTINE  ADAPTED  FRU  4  LINEAR  PR JSRAKMI NG  CODE  LPH-1/  WRITTEN 
OY  .1.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 
***uF.  SCRIPT  ION  nr  PARAMETERS*** 

IP  AH  =  PAS  AMPTER  INDICATING  WHICU  ETA- VECTORS  MATRIX  (ALL  E  OR 
JUST  U  OF  LU  DECUMP)  IS  USED  TO  UPDATE  COLUMN  Y  (INPUT) 

L' 

IMPLICIT  RE AL*4  (  A/C/  F-H/O/ P,  R-W/7.) /  RE AL*8  (8/D/X/Y)/ 

I  INTSGKKM  (T-N/Q) 

INTEGER  JH/KINDAS/LA/Lh/IA/IE 
DOUBLF  PRECISION  E(2000) 

W^AL  A(IOOO) 

C 

COMMON  XLB(122)/XUU(122)/DF/OP,B(60)/X(60)/Y(60)/YTEMP(60)/A/ 

1  E/MSTAT/IObJ/IROWP/ ITCNT/ I MVFRQ/ ITRFRQ/ JCOLP/ 

2  NROW/VC1L/NELEH/NETA/NLELEM/NLETA/NUELEM/NUETA/ JH(60)/ 

3  <INB AS( 122 )/L A (122)/LE(502)/ I A( 1000)/IE(2000) 

c 

HLE  =  HF.TA 
HFE  =  1 

IP  (IPAR.KQ.2)  NFE  =  NLLTA  ♦  1 
IF  (NFE  . GT.  NLf)  RETURN 
DO  1000  IK  =  FFE/NLE 
Ll.  =  LF.(IK) 

KK  =  LP ( I K+l )  -  1 
IPIV  =  1E(LL) 

UY  =  Y(IPIV) 

DY  =  DY/E(LL) 

Y(IPIV)  =  DY 

IF  (KK  .LE.  LL)  GO  TO  1000 
LL  s  ll  ♦  1 
DO  *500  J  =  LL /KK 
IV  =  Tfc(J) 

V(IP)  *  Y(IR)  -  E(J>  *  DY 
bOC  CONTINUE 
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1000  CONTINUE 

RETURN 
FND 

C - - — - - — - - — - — - - - - - — 

SUBROUTINE  BTR  AN 
C 

C  PERFORM  BACKWARD  TRANSFORMATION  ON  COLUMN  STORED  IN  VECTOR  Y 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 

C  b Y  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C 

IMPLICIT  REALM  (A/C/E-H/O/P/R-W/Z)/  REALM  (B/D/X/Y)/ 

1  INTEUERM  (I-N/Q) 

INTEGER  JH/KIM13AS/LA/LE/IA/IE 
DOUBLE  PRECISION  E(2000) 

REAL  A(100C) 

C 

COMMON  XLB(122)/XUB(122)/DE/DP/B(60)/X(60)/Y(60)/YTRMP(60)/A/ 

1  E/MSTAT/ I OBJ/ 1  ROW P/ ITCNT/INVFRQ/ITRFRQ/JCOLP/ 

?  NR3W,  S'  COL  /  NELEl./  NET  A/  NLELEM/  NLET  A/  NUELEM/  NUETA/ JH(60 )/ 

J  KINnAS(122)/L»<122)/LK(502)/IA(1000)/IE(2000) 

C 

1 p  (NETA  .LE.  0)  RETURN 


DO  1000  I  = 

1/NETA 

TK  = 

NETA  -  I  ♦ 

1 

LL  = 

LF(IK) 

KK  = 

LE(IK*1)  - 

1 

IPIV 

=  IF (LL) 

DP  - 

L(LL) 

DY  = 

Y(IPIV) 

DSUM 

=  0. 

I F  (*K  .LE.  LL) 

GO 

LL  = 

LL  ♦  1 

DO  50C  J  =  LL/KK 

IH  =  IE(J) 

D*  =  E ( J ) 

OPROO  =  DE  *  Y( I R) 

DSUM  =  DSUM  ♦  UPROD 
‘jOC  CONTINUE 

oOC  Y(TPIV)  =  ( Li Y  -  OSUM )  /  DP 

10(0  CONTINUE 

RETURN 
E  VD 

C - 

S'ltt-OUTINE  FORNC 

C  FCRM  OBJECTIVE  FUNCTIUN  VECTOR;  IF  BASIS  IS  INFEASIBLE/  SET 

C  OBJECTIVE  FUNCTION  TO  BE  TVFFASimiTY  FORM  FOR  PHASE  I- 

C  C At.Lr-.  FROM  SUBROUTINES  NORMAL  AND  POUNDR. 

r  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 

C  bY  J.  A.  TOMLIN  ( GPFR ATTlJNS  RESEARCH/  STANFORD  UNIVERSITY) 

IMPLICIT  RESLM  <A/C/F-J/C/P/R-R/Z),  RE  AL*8  (B/D/X/Y)/ 

1  P’TEJKkM  (I-v,Q) 

INTEGER  JH/^ISbAS/LA/LE/IA/IE 
DOUHL*  P«  EC  IS  ION  E(200C) 

ALAI.  A(lCOu) 

( 

C  1MMON/C jNSTS/  ZTULZE/ ZTOLPV/ ZTCOST, ZTULSM/ NEGI NF /NEMAX/NRMAX/ QBL/ 
1  NI.F S/ NT T/QA/UT/QF/QN/QSUn,  m/OC/QE/QH/QL/QQ/QR/OM/QG 
COMMON  XLH(122)/XUB(122)/DE/DP, n(EC)/X(60)/Y(60)/YTEMP(60)/A/ 


n  n  n  o  n  n 
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1  F/MSTAT/  IOBJ/  I  RO*P/  1'iCNT/  I  NVFPG/  ITRFRQ/ JCOLP/ 

2  NROW/NCIL/NELtM/UFT  4/  NLFLEM/.'JLET  A/NUELEM/NUET  A/ 

3  RIN3AS(122)/LA(122)/LE(502)/IA(1000)/ IE(2000) 


JH(60)/ 


MSTAT  =  JF 
Y(inBJ)  =  0. 

Dll  30  1  =  1/  NRQW 

IF  ( £  . FQ.  IQBJ)  GO  TO  30 
[CriL  =  JH(  I ) 

IF  (X(I)  .LE.  (XLfl(ICUL)  -  ZTOLZE))  CO  TO  10 
IF  (Y(I)  .  GE.  (XUD(ICOL)  ♦  ZTOLZE) )  GO  TO  20 
Y(I)  =  0. 

Gil  TO  30 


A.  0 

Y(T)  =  i. 

MSTAT 

=  01 

GO  TU 

30 

20 

Y(  I  )  =  -1. 

MSTAT 

=  01 

30 

CONTINUE 

IF  (MSTAT. EG. OF)  Y(IOBJ)  =  1. 
RETURN 

END 


SUBROUTINE  PRICE 

PRICE  OUT  NONBASIC  COLUMNS /  CHOOSE  PIVOT  COLUMN  JCOLP  FOR 
CURRENT  PRIMAL  SIMPLEX  ITERATION.  JCOLP=0  ==>  DUAL  FEASIBLE. 

SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 
BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

IMPLICIT  REALM  (A/C/E-H/O/P/R-W/Z)/  REALM  (B/D/X/Y), 

1  INTPCERM  (I-N/Q) 

INTEGER  JU/KINBAS/LA/Lt/IA/IE 
I MTECER  TPARl/ INCUMB/ IVBNP/ IYTD/ IOBND 
DOUBLE  PRECISION  E(2000) 

REAL  A(1000) 

COMMUV/CUNSTS/  ZTOLZE/ZTOLPV/ZTCOST/ ZTOLSM/ NEGI NF/NEHAX/ NRMAX/QBL/ 
l  NLES/VTMAX/GA/QI/OF/QN/QSUB/QB/UC/QE/QH/QL/OO/UR/QM/QG 
COMMON  XLB(122)/XUB(122)/DE/DP/B(60),X(60)/Y(60)/YTEMP(60)/A/ 

1  E/MSTAT/  IOBJ/  IRUwP/  TTCNT/  INVFRQ/  ITRFRQ/ JCOLP/ 

2  NROW/MCOL/NELEM/NETA/NLELEM/NL£TA/NUELEM/NUETA/JH(60), 

3  KINBAS(122)/LA(122)/LE(502)/IA(1000)/IE(2000) 

CMIN  =  1.E10 
CMAX  =  -1.E10 
DO  1000  J=l/NCOL 

IF  (KINBAS(J)  .GT.  0)  GO  TO  1000 
IF  ( (XUB( J)  -  XLR(J))  .LT.  ZTOLZE)  GO  TO  1000 
CALCULATE  DPR T CF  =  PRICE  OF  BRINGING  COLUMN  J  INTO  THE  BASIS 
DPRICE  =  0. 

LL  =  LA(J) 

KK  =  L A( J+l )  -  1 
DO  500  I  =  LL/KK 
IR  =  I A( I ) 

DE  =  ACT) 

DPRICK  =  DPRICE  ♦  (DE  *  Y(IR)> 

IF  (KINDAS(J)  .EQ.  -1)  GO  TO  600 
IF  (DPRICE  .GE.  CMIN)  GO  TO  1000 
CMIN  =  OPR  ICE 


F 
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JCOLl  =  J 
GO  TO  1000 

600  IF  (DPRICE  .Lfc.  CMAX)  GO  TO  1000 

CVAX  =  DPRICE 
JCOL2  =  J 
1000  CONTINUE 

C 

C  CHOOSE  PIVOT  COLUMN  JCOLP  BASED  ON  PRICES 


c 

IF  (CM IN  .LE. 

-ZTCOST) 

GO  TO 

1500 

IF  (CMAX  .GE. 

JCOLP  =  0 
RETURN 

ZTCOST) 

GO  TO 

2000 

1500 

IF  (CMAX  • GE. 

ZTCOST) 

GO  TO 

2500 

1600 

JCOLP  =  JCOL1 
RETURN 

/000 

JCOLP  =  JCOL2 
R-TURN 

2500 

IF  (A"’3(CMIN) 

-  CMAX) 

2000,2000,1600 

. 

END 

SUBROUTINE  CilUZK 

C 

c  PERFORM  FIN-RATIO  TEST  FOR  PIVOT  COLUMN  JCOLP  DETERMINED  IN 

C  SUBROUTINE  PRICF,  THEN  SELECT  PIVOT  ROW  IROWP  FOR  CURRENT 

C  PRIMAL  SIMPLEX  ITERATION. 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1,  WRITTEN 

C  BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH,  STANFORD  UNIVERSITY) 

C 

IMPLICIT  REILM  ( A,C,E-H,0,P,  R-W,Z),  RE  AL*8  (B,D,X,Y), 

1  IHTECKSM  ( I-N,Q) 

INTEGER  J!I,KI'1B4S,LA,LE,IA,IE 
INTEGER  IPAPT, INCUMB, IVBND, IVID, IOBND 
DOUBLE  PPECTSTUN  E(2000) 
kEAL  A(1C00) 

C 

cn  'MON/crNSTS/  ZTOLZE, ZTOLPV, ZTCOST, ZTOLSM/ NEC I NF , NEMAX, NRM AX, QBL, 
l  NL;<S,N?MAX,QA,QI,QF,QN,QSUB,QB,QC,WE,QH,QL,QO,WR,QII,QC 
CnvvnN/R.'LISl/  DFPART(60),RF.VBND,T.ICV4L,  I  COL,  IVAL,IDIR,  IPART(122), 
1  T  VCMP (  130),  rVbND(500),  IV ID (500),  IOBND ( 500 ),NPIVOT,  IPTYPE/IFEAS 
C  WON  XLH(122),XUB(122),DE,DP,B(60),X(60),Y(60 ), V TEMPI 60 ), A, 

1  E  ,  *’ STAX,  rr'LJ,IFOWP,lTCNT,INVFRQ,ITRFKQ, JCOLP, 

?  Nt-’i"/',  NOOL,  \FLrM,  \ETA,  NLELEM,  NLET  A,  NUELEM/NUETA,  .JH(60  ), 

3  K  T  VMS  (122), L 4(122 ),LE (5  02),  I  »  (  1000),  IE (2000) 

(’ 

r  V  'CriR  V  CONTAINS  FORWARD  TRANSFORM  0*  COLUMN  JCOLP. 

C  JP  =  WIN.  RA1IU  SO  FAR.  IT  IS  PASSED  TO  UPBFTA  FOR  THE  PIVOT  STEP. 

C  also  PASSED  TO  UPBLTA  IS  DE  C  =  X(IPOWP)  AFTER  PIVOT  STEP  3. 

C  \n  I  VOT  =  0  I Fr  JCOLP  IS  NON? AS IC  (AT  OPPOSITE  BOUND)  AFTER  PIVOT. 

C. 

re1  (F INR  AS  (JCOLP )  .EQ.  -1)  SO  TO  100C 
C 

C  I'TOPHS  V  ARI 4PLr  AT  LOWER  BOUND?  COMPUTE  MIN  RATIO  DP 

C 

D r>  =  1.H0 
TO  500  1  =  1,  NMD* 

If  (I  .FQ.  I DBJ)  GD  TO  500 
ICOL  =  JP(!) 

it  ( Y ( I )  .ST.  ZTOLPV)  GO  TO  100 
IF  ( Y ( I )  .LT.  -ZTOLPV)  CID  TC  200 
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GO  TO  500 
C 

IOC 


c 

200 


250 

500 


600 
C 

C  INCOMING  VARIABLE  AT  UPPER  BOUND?  COMPUTE  MAX  RATIO  DP 
C 

1000  DP  =  -1.E10 

DO  1500  I  =  l/I.ROW 

IE  (I  .EG.  IOBJ)  GO  TJ  1500 
ICUL  =  JH( I ) 

IF  (  Y  (  T  )  .GT.  ZTOLPV)  GO  TO  1100 

IF  (V(l)  .LT.  -ZTOLPV)  GO  TO  1200  j 

GO  TO  1500  / 

C  POSITIVE  COEFFICIENT  A(I,JCULP) 

110C  IF  (X(I)  .GT.  (XUB(ICUL)  «■  ZTOLZE))  C-0  TO  1500 

OF  =  (XCI)  -  XUB(rCUL))/Y(I) 

IF  (DF.  .LE.  DP)  GO  TO  1500 
IPTYPF  =  -1 
GO  TO  1250 

C  NEGATIVE  COEFFICIENT  A(I,JCOLP) 

120C  IF  <X(I)  .LT.  (XLB(  ICOL)  -  ZTOLZE) )  GO  TO  1500 

PE  =  (X(I)  -  XLb( ICOL)  )  / Y ( T ) 

IF  (DE  .LE.  DP)  GO  TO  1500 
IPTYPE  =  0 
1250  DP  =  CE 

IROKP  =  T 
1500  CONTINUE 

DF  =  DP  ♦  XUP(JCOLP) 

IF  (DE  .GT.  XLE(JCOLP))  GO  TO  1600 
DP  —  XLB(JCULP)  -  XUB( JCllLP) 

NPIVOT  =  0 
RETURN 

1600  NPIVOT  *  1 
RETURN 
END 

SUBROUTINE  UP  ETA 
C 

C  FORM  NEh  ETA-VECTORS  FOR  PRWUCT  FORM  OF  BASIS  INVERSE 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 

C  BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 


POSITIVE  COEFFICIENT  A(l/JCOLP) 

IF  (X(I)  .LT.  (XLB(ICOL)  -  ZTOLZE))  GO  TO  500 
Pt  =  (X(I)  -  XLF»(ICJL))/V(T) 

Tr  (DE  .Gh.  DP)  GO  TO  500 
IPTYPE  =  0 
GO  TO  250 

NEGATIVE  COEFFICIENT  A(I,JCJLP) 

IF  (XCI)  .GT.  (XUR(ICOL)  ♦  ZTOLZE) )  GO  TO  500 
OF  =  (X(I)  -  XUb( ICOL) )/ Y( I) 

IF  (DF.  .Gc.  DP)  GO  TO  500 
IPTYPE  =  -1 
np  =  OF 

IPO¥P  =  I 
CONTINUE 

D  •  =  DP  ♦  XLb(JCULP) 

IF  (DE  .LT.  XUP(JCOLP))  GO  TO  600 
DP  =  XUB(JCOLP)  -  XLB(JCULP) 

NPIVOT  =  0 
RrTURN 
NPIVOT  =  1 
RETURN 


j 


•">  a  n  cj  n  n 
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C 

IMPLICIT  RE AL*4  (  k,C,E-H,  0,  ?,  R-W,  Z  ),  REAL  *3  (B/D/X/Y)/ 

1  INTEGERM  (I-N/Q) 

INTEGER  JH/K1NBAS/LA/LE/IA/IE 
DOUBLE  PRECISION  E(2000) 

REAL  A(1000) 

C 

COMMON/CCNSTS/  ZTOLZE/ZTOLPV/ZTCOST/ZTOLSM/NEGINF /NEMAX/NRMAX/QBL/ 
1  NiLES/NTMAX/QA  /  QI/QF/QN/  QSUB/QB,  QC/QE/  QH/QL/ QO/QR/QX/QG 
COMMON  XLB(122)/ XUB( 122)/DE/DP/B(60)/X(60)/Y(60)/ YTEMP(60), A, 

1  F/MSTAT/  IOBJ/1KOWP/  ITCNT/INVFPa/ITRFRQ/JCOLP/ 

2  NRON/NCOL/NFLEM/NETA/NLFLEM/YLSTA/NtJELEM/NUETA/  JH(60), 

1  K INB AS(1 22)/LA(l 22)sLE(502)/ I A( 1000)/ T£( 2000) 

C 

NELEH  s  NELEM  ♦  1 
IE(NELEM)  -  1R0WP 
E(NE'LEM)  =  Y(IROWP) 

DO  1000  I  =  l/.NROW 

IF  (I  .EQ.  IROWP)  GO  TO  1000 

IF  (OABS(Y(I))  .LE.  ZTOLZE)  GO  TO  1000 

NELE4  =  NELEM  ♦  1 

IE(NELEM)  =  T 

E(YELEM)  =  Y( I ) 

10C0  CONTINUE 

N  K  T  A  =  NET A  ♦  1 
LE(NETAd)  =  NELEM  ♦  1 
RETURN 
EMU 

C - 

SUBROUTINE  SHIFTKdQLG/  INEW) 

REARRANGE  DATA  STORAGE;  USED  BY  SUBROUTINE  INVERT 
SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1#  WRITTEN 
BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 
••♦DESCRIPTION  OF  PARAMETERS*** 

I'JLD/INEM  =  PARAMETERS  INDEXING  STORAGE  LOCATIONS  IN  WHICH 
C  OATA  IS  TO  BE  TRANSFERRED  (INPUT) 

v.. 

IMPLICIT  »EAL*4  ( A/ C/E-H/O/P/ R-W/Z) /  RE AL*8  (B/D/X/Y)/ 

1  INTEGER*4  (I-N/Q) 

INTEGER  JH/XINBAS/LA/LE/IA/IE 
DOUBLE  PRECISION  E(2000) 

REAL  A(1C00) 

C 

C  IMMO’J/CJNSTS/  ZTOLZE/ ZTULPV /  ZTCOST/  ZTOLSM/NEGINF/  NEMAX/NRN  AX/4BL/ 
1  NLFS/NTMAY/OA/ai/QF/QN/QSUB/QR/QC/ilE/OH/QL/OO/OR/aM/QG 
COMMON  XLR(12?)/XUB(122)/DS/DP/R(60)/Y(60)/Y(60)/YTEMP(60)/A/ 

1  E/EST AT/ IObJ/ IKQWP/ ITCNT/ l YVFPQ/ IT RFRQ/ JCOLP/ 

>  NRjW/NCnL/*  EcFM/NET A/  NLELKM/NLET A/NUELEM/NUET A/ JH(60.)/ 

<  KrNBAS(12?)/LA(122)/LE(50?)/I »( 1 000)/ IE( 2000) 

C 

DIMENSION  B  ARRAY (240) 

EQUIVALENCE  (PARRAYd  )/P(l  )) 

IFO  =  (IDLE  -  1)  *  NRMAY 
IEM  =  (ISEk  -  1)  *  NRMAX 
DO  1000  i  =  1/NRUW 

OARRAY(lEN  ♦  I)  =  BARR  AY  (IFO  ♦  I) 
continue 

RETURN 

r  no 


1000 


n  n  ~i  nnocro 
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SUBROUTINE  INVERT 

C 

C  COMPUTE  INVERSE  OF  CURRENT  BASIS  BY  LU  DECOMPOSITION 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPH-1/  WRITTEN 

C  MY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C 

IMPLICIT  REALM  (A/C/F-H/Cl/P/P-W/7.)/  REALM  (R/D/X/Y)/ 
l  TKTFUERM  (I-N/Q) 

INTEGER  Jrt/KINBAS/LA/LF/IA/ll' 

DOUBLE  PRECISION  F(2C00) 

REAL  A(1000) 

C 

COmKON/CONSTS/  ZTOLZE/ETOLPV/ ZTCOST/ZTOLSM/NEGINF/NEM AX/NRMAX/QBL/ 
l  MLES/NTMAX/QA/UI/QF/QN/QSIlU/QB/QC/QE/QH/UL/QO/aR/QM/UG 
COMMON  XL8(122)/XU8(122)/DE/DP/B(60)/X(60)/Y(60)/YTEMP(60)/A/ 

1  F/VSTAT/ 1DPJ/ IROMP, ITCNT/ I NVFRQ/ ITRFRQ/JCOLP/ 

2  NROW/ NCOL/ NKLEM/NETA/NLELEM/NLETA/NUELEM/NUETA/ JH(60)/ 

3  KTNBAS(122)/LA(I22)/LE(502)/IA(1000)/ IE (2000) 

Cu MMON/T I  MENS/ ITOT/T STORE/ TIMFLP/TIMEOR/TIMEOC/TIMIWV 
C 

INTEGER  MREG/HREG/ VRKG 
0 IMFNS  IQ  S'  MRLG  (60  )/HRKG(oO  )  /  VREG(60) 

EQUIVALENCE  (MRFG(1)/YT£MP(1) )/ (HREG<  1 )/ YTEMPOl ) ),(VREG(1),X<1)) 

I TIME  =  IHPTIM(l) 

SET  PARAMETERS 

NET  A  =  0 
NLETA  =  0 
NUETA  =  0 
NELEM  =  0 
NLELEM  =  0 
NUELEM  =  0 
N ABOVE  *  0 
LE ( 1)  =  1 
Kll  =  C 

LH4  =  NROW  ♦  1 

PUT  SLACKS  AND  ARTIFICIALS  IN  PART  4  AND  REST  IN  PART  l 

00  100  I  =  l/NRP* 

IF  ( JH( I )  .GT.  NROW)  GO  TO  50 
Ln4  =  LF4  -  1 
MMEGCLR4)  =  JH( I ) 

VREC(LR4)  =  JH(I) 

GO  TO  90 

LO  KR1  *  FR1  ♦  1 

VREG(KRl)  =  JH( 1 ) 

90  HREG(I)  =  -1 

JH( I )  =  0 
100  CONTINUE 

FR3  =  LR4  -1 

DO3200LI  *  LR4/NROW 
IF  *  SREC(I) 

HREC(IR)*  0 
JH(IP)  =  IR 
KIMBAS(IR)  *  IR 
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200  CONTINUE 

C 

C  PULL  OUT  VECTORS  BELOW  BUMP  AND  GET  ROk  COUNTS 
C 

NBNONZ  =  NfiOW  -  LR4  ♦  1 
IF  (KR1  .EG.  0)  GO  TO  1190 
J  =  1 

210  IV  =  VHEG(J) 

LL  =  LA(IV) 

KK  =  LA( I V+l )  -1 
IRCNT  =  0 
DO  220  I  =  LL,KK 

NriNONZ  =  NBNONZ  ♦  1 
IP  =  IA( T ) 

IF  (HRFG(IR)  .GE.  0)  GO  TO  220 
TKCNT  =  IRC NT  ♦  1 
HKEG(IR)  =  HREG(IR)  -  1 
IRP  =  IR 

220  CONTINUE 

IF  (IRCNT  -  1)  230/ 250/ 300 
230  *  ’ITE(21/8000) 

9000  F-|R>IAT(16HKATRIX  SINGULAR  ) 

KThBAS(IV)  =  0 
VNEC(J)  =  VPEG(KRl) 

KR1  =  KR1  -  1 
IF  (J  .Gf.  Kkl )  GO  TO  310 
G1  TO  210 
C 

250  VPEG(J)  =  VREG(Kkl) 

KR1  =  KR1  -  1 

LR3  =  LR3  -  1 

V>EG(LR3>  =  IV 

KRFG(LW3)  =  IRP 

HRES(JRP)  =  0 

JM(IRP)  -  IV 

KTMBAS ( IV )  =  IRP 

IF  (J  .  GT.  KR1)  GO  TO  310 

r,l  TO  210 

300  IF  (J  . GrJ .  KR1)  GO  TO  310 
J  =  J+l 
G  I  TO  210 
C 

C  f,,JLL  UUT  REMAINING  VECTORS  ABOVE  AND  BELOW  THE 
C  iJUUP  AND  FSTAILISH  MERTT  COUNTS  OF  COLUMNS 
r 

310  KVREM  =  0 

IF ( KP 1  • fcQ.  0)  GO  TO  1190 
J  =  1 

320  IV  =  VPSG(J) 

LL  =  LA(TV) 

KK  =  LA(IV4l)  -  1 
IRCNT  =  0 
L'T  BOO  I  =  LL/KK 
IR  =  I A( I ) 

IF(H»LG(IK)  . NE.  -2)  GO  TO  400 

C 

c  pivot  a »uv f;  bump  (part  of  l) 

c 

'lAHOVr  v  MiJOVF  ♦  1 
I  h  Q  *  P  =  I  U 
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CALL  UNPACK(IV) 

CALL  WRaTA 
HL’-:TA  =  NETA 
JH(IR)  =  I? 

KINBAS(IV)  =  IR 
VREG(J)  =  VREG(KRl) 

KR1  =  KR1  -  1 
NVREM  =  NVREM  ♦  1 
HREG(IR)  =  IV 
GO  TO  940 
C 

400  IF  (MREG(IR)  .GE.  0)  GO  TO  900 

IRCNT  =  IRCNT  ♦  1 
IRP  =  IR 

bOC  CONTINUE 

C 

IF  (IRCNT  -  1)  810/900/1000 
810  WRITE (21/ 8000  ) 

KINflAS(IV)  =  0 
VREG(J)  =  VR£G<KR1) 

NVREM  =  NVREM  ♦  1 
KR1  =  KR1  -  1 

IF  (J  .GT.  Kkl)  GO  TO  1010 
GT  T'J  320 
C 

C  PUT  VECTOR  BELOW  BUMP 

r 

900  VREG(J)  =  VREG(XRl) 

NVH&M  =  NVRFh  ♦  1 
KR1  =  KR1  -  1 
LR3  =  LRJ  -  1 
VRKG(LK3)  =  IV 
MREG(Lk3)  =  IRP 
HREG(IRP)  =  0 
J'K  IRP)  =  IV 
KTMRAS(IV)  s  IRP 
C 

C  CHANGE  ROW  COUNTS 
C 

940  00  950  11  =  LL/KK 

IIR  =  I A ( 1 1 ) 

IF  (HPEG(IIR)  .GE.  0)  GO  TO  950 
HREG(IIR)  =  HRFG(IIR)  ♦  1 


950 

CONTINUE 

IF 

(J  .GT. 

m 

) 

GO 

TO 

1010 

CO 

TO  320 

100^ 

IF 

(J  .GE. 

rfcl 

) 

GO 

TO 

1010 

J  : 

=  J>1 

cn 

TO  320 

101C 

IF(NVREM  . 

GT. 

0) 

GO 

i  TO 

310 

C 

C  GET  MERIT  COUNTS 
C 

1020  IF  (KR1  .EQ.  0)  GO  TO  1190 
DO  1100  J  =  1/KR1 

IV  =■  VPEG(J) 

LL  =  LA(TV) 

Kf  r  LA( I V+l )  -  1 
IMCNT  =  0 
DO  1050  I  =  LL/KK 
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IR  =  IA( I ) 

IF  (HRkC(IR)  .GE.  0)  GO  TO  1050 
IMCNT  =  IMCNT  -  (HREG(IR)  *1) 

1050  CONTINUE 

MREG(J)  =  IMCNT 

1100  CONTINUE 

C 

C  SORT  COLUMNS  INTO  MERIT  ORDER  USING  SHELL  SORT 
ISO  =  1 

1106  IF  (KR1  .LT.  2*ISD)  GO  TO  1108 
ISO  a  2*1  SI) 

GO  TO  1106 
1106  ISO  a  ISO  -  1 
C  F MU  OF  INITIALIZATION 

1101  IF  (ISO  .L£.  0)  GO  TO  1107 
I3K  =  1 

1102  ISJ  a  rsK 

ISL  a  ISK  ♦  ISO 
ISY  a  MPEG(ISL) 

t  C7  -  VDh'(T  MCI  \ 

1103  IF  ( IS Y  .LT.  MkEG(ISJ))  GO  TO  1104 
1U6  ISL  =  ISJ  ♦  ISD 

>1  s< c G (  I SL)  =  ISY 
VUEG(ISL)  =  1SZ 
ISK  =  ISK  ♦  1 

IF  ((TJK  ♦  ISO)  .LE .  KR1 )  GO  TO  1102 
ISO  a  (ISO  -  1)  /  2 

cn  rn  noi 

11C*  I^L  =  ISJ  ♦  I  SO 

M’cG(IjL)  a  HRLG(IS.J) 

VPfcG(TSL)  a  W^tGdSJ) 

ISJ  =  ISJ  -  ISO 
IF  (ISJ  .GT.  0)  GO  TO  1103 
GO  TO  11 1)5 
1 1 C  7  CONTINUE 

r 

C  uVLi  l)F  SOKT  ROUTINE 
C  P'JT  ">UT  a  BLOW  BUMP  ETAS  (PART  OF  U) 

C 

119C  NSLCK  =  0 
INFLOW  =  0 
N EL  AST  =  NL»'AX 
«TL AST  a  NTVAX 
LF(NTLA  >T  ♦  1)  =  NEL AST  ♦>  1 
Lk  =  LW  3 

IF  (LF3  .Gk.  L*4>  LR  a  LR4 
IF  (Lk  .IT.  NRO*>  GO  TO  2050 
J<  =  VHOw  ♦  I 
on  2000  jja  LP/NROR 
JK  =  JK  -  1 
IV  a  VPF.G(JK) 

T  a  MkEG(JK) 

NHELOk  =  NBKLOM  ♦  l 
IF  (IV  .GT.  NRUk)  GO  TO  1200 
SSLC*  =  NSLC*  ♦  1 
1200  LL  =  LA (TV) 

K ft  a  LA( T  ¥*1 )  -1 
IF  (KK  .GT.  LL)  GO  T1  1J00 
IF  ( A'iS(  l(LL)  -  1.)  .LE.  ZTOLZE)  GO  TO  2000 


1250 
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1 3 C 0  N'M-'TA  =  NUFTA  ♦  1 

DU  1400  J  =  LL/KK 
IR  =  I  A(  J) 

IF  (IR  .EQ.  I ) GO  TO  1390 
lFUtlLAST)  =  I.n 
E ( NEL AST )  =  A(J) 

NEL AST  =  MELAST  -  1 
MJELEM  =  NUELFM  ♦  1 
GO  TO  1400 

1390  KP  =  KJ) 

14C0  CONTINUE 

l'(NFLAST)  =  £ 

E  (  Ni,L  AST  )  =  £P 
LE(NTLAST)  =  NELAST 
'isLAFT  =  KELAST  -  1 

•JTLAST  =  NTL AST  -  1 

NUELFM  =  NUtLE-1  ♦  1 

20C  0  CONTINUE 

2050  IF(K*1  .  EQ.  0)  GO  TO  3500 
C 

C  Dfl  L-U  DECOMPOSITION  OF  BUMP 
C 

DO  3000  J  =  1,KP1 

IV  =  VRfcG(J) 

CALL  UNPACK(IV) 

CALL  FTKAN(2) 

IROWP  =  0 
I  KCMTN  =  -999999 
DO  2100  I  =  1 / NROW 

IF  (DABS(Y( I) )  .LE.  ZTOLPV)  GO  TO  2100 
IF  (HREG(I)  .GE.O)  GO  TO  2100 
IF  (HREG(I)  .LE.  IRCMIN)  GO  TO  2100 
IK CM IN  =  HKEG(I) 

IROWP  =  I 

2100  CONTINUE 

TF  ( IROWP  .GT.  0)  GO  TO  2150 
WRITF<21,8000) 

KINBAS(IV)  =  0 
GJ  TO  3000 

2150  INCH  =  HREC( IROWP)  ♦  3 

C 

C  WRITE  L  AND  U  ETAS 
C 

IF  (J  .EQ.  KR1)  GO  TO  2160 
NKLEM  =  NELEM  ♦  1 
IK( NELEM)  =  IROWP 
F( NELEM)  =  Y( IROWP) 

2160  DO  23C0  I  =  1/NROW 

IF  (I  .EQ.  IROWP)  GO  TO  2300 
IF<DABS<YU))  .LE.  ZTOLZE)  GO  TO  2300 
IF  (HREG(I)  .GE.  0)  GO  TO  2200 
C 

C  f.  ETA  ELEMENTS 
C 

NELEM  *  NELFM  ♦  1 
lE(NELEM)  =  I 
t(  NELEM)  =  YU) 

GO  TO  2300 
C 

C  U  ETA  ELEMENTS 
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C 

2200  lE(NELAST)  =  I 

E( NELAST)  *  Y(I) 

N EL  AST  =  MEL  AST  -  1 

Ml  EL EM  =  NUELEM  ♦  1 

2300  CONTINUE 

C 

JH(IRONP)  =  IV 
KINBAS(IV)  =  IRONP 
NUETA  =  NUETA  ♦  1 
IE(NELAST)  =  IRQMP 
IF  <J  .ME.  KR1 )  GO  TO  2330 
E(NELAST)  =  Y(IROWP) 

GO  TO  2340 
2330  E(HELAST)  =  1. 

NETA  =  NET A  ♦  1 

LE(NETA+1)  =  NELEK  ♦  1 
2340  NUELEM  =  NUELEM  ♦  1 

LE(NTLAST)  =  NELAST 
NELAST  =  NELAST  -  1 
NTL  AST  =  NTL AST  -  1 
C 

C  UPDATE  ROW  COUNTS 
C 

DO  2350  I  =  1,NR0N 

IF  <D ABS(Y( I) )  .LE.  2T0LZE)  GO  TO  2350 
IF  (HREG(I)  .GE.  0)  GO  TO  2350 
HR EG ( I )  =  HREG(I)  -  INCR 
IP  (HREG(I)  .GE.  0)  HREG(I)  *  -1 
2350  CONTINUE 

HREG(IPCKP)  =  0 
3000  CONTINUE 

C 

C  MERGE  L  AND  U  ETAS 
C 

3500  NLETA  =  NETA 

NETA  =  NLETA  ♦  NUETA 
NLELEM  =  NLLEM 
NFLEM  =  NLELEM  ♦  NUELEM 
IF  (NUE.LEM  . LO.  0)  GO  TO  3550 
CALL  SHFTE 
C 

C  INSERT  SLACKS  FOR  DELETED  COLUMNS 

C 

3550  LJ  3500  I  =  1,NR0M 

IF  (  J  E  (  I )  .NF..  0)  GO  TO  3600 
JHU)  =  I 
1*0*?  =  I 

CALL  UNPACK! IROWP) 

CALL  FTRAN(l> 

CALL  ARFTA 
3600  CONTINUE 

C 

C  UPDATE  X 
C 

CALL  URDATA 

L  I  TIME  2  s  IHPTIN(l) 

u  TIMINV  =  TIMINV  ♦  (I T! MF2- 1  TIME )/ 1 00000. 

RETURN 
(•  NO 
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SUBROUTINE  UNPACK(IV) 

C 

C  EXPAND  COMPRESSED  MATRIX  COLUMN  AND  STORE  IN  VECTOR  Y 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-I,  WRITTEN 

C  PY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C  *** DESCRIPTION  OF  PARAMETERS*** 

C  IV  =  PARAMETER  INDEXING  COLUMN  TO  BE  EXPANDED  (INPUT) 

C 

IMPLICIT  REALM  ( A/C/F-H/O/P/R-W/Z)/  R£AL*8  (B/D/X/Y)/ 
t  INTEGER**  (I-N/Q) 

INTEGER  JH/KINBAS/LA/LE/IA/IE 
DOUBLE  PRECISION  E(2000) 

RIAL  A(1000) 

r 

COMMON  XLd(122)/X'J'i(122)/ni;/DP/B(60)/X(60)/Y(60)/YTEMP(60)/A/ 

1  F/MSTAT/  IOBJ/  IRONP/  ITCNT/  I NVFRQ/  ITRFRQ/ JCOLP/ 

2  NROW/NCOL/NELEM/NETA/NLELEM,NLETA/NnELEM/ NUE T A/ JH( 60) / 

3  KINBAS(12?)/LA(122)/LF.(502)/IA(1000)/I£(2000) 

C 

DJ  100  I  =  1/NROW 
Y(I)  =  0. 

100  CONTINUE 

LL  =  L  A ( l V ) 

KK  =  LA( IY+1)  -  1 
DO  200  I  =  LL/Kff 
IP  =  I A( I ) 

Y(IR)  =  A ( I ) 

200  CONTINUF 

RETURN 
END 

C - 

SUBROUTINE  SHFTE 
C 

C  SUBROUTINE  FOR  INVERT 

C  SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 

C  BY  J.  A.  TOMLIN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C 

IMPLICIT  RE ALM  ( A/C/E-H/O/P/R-W/Z)/  REAL*8  (B/D/X/Y)/ 

1  INTEGEKM  (I-N/Q) 

INTEGER  Jli/FINdAS/LA/LE/IA/IF. 

DOUBLE  PRECISION  F.(2000) 

REAL  A(IOOO) 

C 

COHMUN/CUNSTS/  ZTOLZE/ZTOLPV/ZTCOST/ZTULSM/NEGI NF/NEMAX/ NRMAX/OBL/ 
1  NLES/NTMAX/aA/QI/UF/UN/QSUH/OB/QC/QE/QH/QL/QO/QR/QM/QG 
COMMON  XLB(122)/XUR(122)/DE/DP/B(60)/X(60)/Y(60)/YTEMP(60)/A/ 

1  E/MSTAT/IOBJ/T ROMP/ ITCNT/ I NVFRQ/ ITRFRQ/ JCOLP/ 

2  NROW/NC0L/NFLFM,UETA/NLELFM/NLETA/NUELEM/NUETA/JH(60)/ 

3  KINPAS( 122)/L A(122)/LE(502)/I A( 1000),IE( 2000) 

C 

C  SHIFT  IE  AND  F  OF  U  ELEMENTS 
C 

NF  =  MEMAX  -  NUELEN  ♦  1 
INCR  a  0 

DO  1000  I  a  NF/NcMAX 

INCR  a  [NCR  4  1 

Ir. (NLELfcM  4  INCR)  a  IE(I) 

E(NLELEM  4  INCR)  a  £(I) 

1000  CQNTINOg 


o  o  r.  n  o  o  n  r>  n  n  m  ►-  n  o  r,  ^  n  r»  o  o 
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C 

IOIF  =  N'MAX  -  NLELEM  -  NUELKM 
NF  =  NTMAX  -  MJFTA  ♦  1 
Tiiru  c 

DO  2000  1  =  NF/NTMAX 

raru  s  rncp  +  1 

LEOILETA  ♦  INCR)  s  L£(T)  -  IDIF 
2000  CONTINUE 

LEI NETA*1 )  s  NELEM  ♦  1 

RETURN 

END 

SUBROUTINE  UpBETA 

UPDATE  RTGHT-H AND  SIDES  TO  REFLECT  NEW  BASIS  RESULTING  FROM 
CURRENT  SIMPLEX  PIVOT 

SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 
BY  J.  A.  TOMLIN  ( OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

IMPLICIT  REALM  ( A/C/E-H/O/P/ R-W/Z)/  RfcAL*8  (B/O/X/Y)/ 

1  INTEGERM  (I-N/Q) 

INTEGER  JH/KINBAS/LA/LE/IA/1E 
INTEGER  I PART/ INCUMB, IVBND/IVID/IOBND 
DBUBLE  PRECISION  E(2000) 

HEAL  A(1000) 

COMMON/BBLIST/  DPP ART< 60 ) , RE V BND, T hCV AL, 1COL, I V AL, IDIR, IP ARTC122 ) , 
1  INCUMB(130)/IYBND(500)/IVID(500)/IOPND(500)/NPIYOT/IPTYPE/IFEAS 
COMMON  XLB(122)/XUB(122)/DE/DP/B(60)/X(60)/Y(60)/YTEMPC60)/A/ 

1  E/RSTAT/IOBJ/IROWP/ITCNT/INYFRQ/ITRFRQ/JCOLP/ 

2  NR0V,!NCnL/NELEM/NFTA/NLELEM/NLETA/NUELEM/NUETA/JH(60)/ 

3  K I NBAS! 122)/LA(122)/LE(502)/I 4(1 000)/ IE ( 2000) 

DO  1000  Isl/NROW 
000  X(I)  =  /(I)  -  Y(I )*DP 

IF  (NPIYOT  .EG.  1)  GO  TO  2000 
Kri"A<5(J<'OLP)  3  -(KINBAS(JCOLP)  1) 

RETURN 

000  A(IPOWP)  =  DE 

IVO^T  s  JH(IMOWP) 

KINBAS(JCOLP)  s  IROWP 
KT NBA  3( IVOUT)  s  IPTYPS 
JU(IROWP)  3  JCOLP 
RETURN 
».Nf> 


SUBROUTINE  NOPmAL(ITSINV) 

THIS  IS  THE  MASTER  PROGRAM  FOR  LINEAR  PROGRAMMING  COMPONENT 
( HEY ISfcD/PK f MAL-SINPLEX  METHOD)  OF  BR ANCH-AND-BOUND  ROUTINE. 
SUBROUTINE  ADAPTED  FROM  LINEAR  PROGRAMMING  CODE  LPM-1/  WRITTEN 
'■V  J.  A.  TOMLIN  (OPERATIJNS  RESEARCH/  STANFORD  UNIVERSITY) 
••^DESCRIPTION  OF  PARAMETERS*** 

1TSINV  3  NUrtBER  OF  SIMPLEX  ITERATIONS  SINCE  LAST  BASIS 
INVERSION  (INPUT/OUTPUT ) 

IMPLICIT  REALM  <  A/C/E-H/O/P,  R-N/Z)/  REALM  (B/D/X/Y)/ 

1  iNTMKSM  (I-N/O) 

INTEGER  JH/MNBAS/LA/LE/IA/IE 
INTEGER  IPART/INCUMH/IVBND/IVID/IOPNO 


UOliBLK  PRECISION  E(2000) 

REAL  A( ICOO) 

C 

CDMMuM/C'-NSTS/  ZTOLZK,ZTOLPV, 7TC0ST, ZTOLSM,  NEGl NF, NEMAX, NRM AX,QBL 
1  NLES,N?MAX,QA,QI,QF,QN,QSU2,QB,QC,QE, QH, QL,  QO,QR,uM,QG 
COMMON/ BKL I  ST/  f)F?ARr(60),REYBNO,INCVAL,  ICOL,  IVAL,IOIR,  IP  ART  (122) 
1  INCUMB(  130  )/  IVBND(500),IYri)(500),I  QBND( 500 ), N PIVOT,  IPTYPE,  I  PEAS 
COMMON  XL3(122),XUB(122)/DE,DO/B(60)/X(60)/Y(60),YTENP{60),A# 

1  R, MSTAT,  IOB  J,  I ROWP,  ITCNT, I NVFRQ, ITRRRQ, JCOLP, 

2  NRO*i,NCOL,NFLEM,  NETA,NL..L2  i,HLETA,NUe;L£M,NUETA,JH<60), 

3  KINBAS(122),LA(122),LE(502),IA(1000),IE(2000) 

r 

17  (I’SIfJV  .LT.  1NVFRQ)  GO  TO  1500 
1000  CALL  INVERT 
ITSINV  =  0 

c 

C  3IMDLFX  CYCLE 
C 

1500  CALL  FOPMC 
CALL  BTRAN 
CALL  PRICE 

IF  (JCOLP  .GT.  0)  GO  TO  3000 
IF  (HSTAT  • EQ.  QI)  GO  TO  2000 
MSTAT  =  1BL 
RETURN 

2000  MSTAT  =  QH 
RETURN 

r 

C  PIVOT  ON  COLUMN  JCOLP. 

C 

30C0  CALL  UNF  ACK( JCOLP) 

CALL  PTR AN (1 ) 

CALL  CWJZR 

CALL  OP 3 ETA 

ITCNT  =  ITCNT  ♦  1 

I  ▼  3 1  'J  V  =  ITSINV  ♦  1 

IF  (NPIVDT  .EQ.  0)  GO  TO  -1010 

IF  (NELEm  .GT.  (NEHAX-NRCJW))  GO  TO  1000 

CALL  WR7  f  A 

1010  IP  (JTSI.W  .GE.  INVFKQ)  GO  TO  1000 
IF  (ITCNT  . GE.  ITRFPQ)  RETURN 
m  TO  1500 
KTU 

S'JbROU  TINE  BANDB(INITBD) 

C 

C  MASTER  PROGRAM  FOR  BRANCH-AND-BOUND  INTEGER  PROGRAMMING 

C  ROUTINE.  ALSO  SERVES  AS  MASTER  PROGRAM  FOR  REOPTIMIZATION 

C  VIA  REVISED  DUAL-SIMPLEX  METHOD  AFTER  A  FORWARD  BRANCH. 

C  SUBROUTINE  ADAPTED  FROM  INTEGER  PROGRAMMING  CODE  BB,  WRITTEN 

C  BY  GARY  A.  KOCHMAN  (OPERATIONS  RESEARCH,  STANFORD  UNIVERSITY) 

C  ***UESCRTPT10N  OF  PARAMETERS*** 

C  INITBQ  =  INITIAL  LOWER  BOUND  ON  MAX.  OBJECTIVE  VALUE  (INPUT) 

C 

IMPLICIT  REAL*!  (A,C,F-H,0,P,R-W,Z),  REAL*8  (R,D,X,Y), 
l  INTEGER*4  (!-M,Q) 

INTEGER  JH,KTN'i»S,LA,LK,IA,IK 
INTEGER  IPART, INCUMB, I V3ND, 1 VID, IOBND 
DOUBLE  PRECISION  F.(2000) 

REAL  A( 1000 ) 
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C 

COHMON/CUNSTS/  ZT0LZB/7.T0LPV/ZTCQST/ZT0LSM/NEGINF /NEMAX/IRMAX/QBL/ 
1  NLES/NTMAX/QA/QI/QP/QN/QSUB/QP/QC/QE/QH/QL/QQ/QR/ QM/QC 
COMMON/BBL IST/  DFPART(60)/RFVBND/ INCV AL/ ICOL/ IVAL/IDIR/ IP ART (122)/ 
l  IMCUMR(130)/IV«NO(500)/IVID(500),IOBND(500)/NPIVOT,IPTYPE/IPEAS 
COMMON  XLB{122)/ XUB(122)/D£/DP/ ft(60)/X(60)/Y(60)/¥TEMP(60)/ A# 

1  K/MSTAT/ IOLJ/IROWP/ITCHT/ TNVFRQ/ ITRFRQ/ JCOLP/ 

2  NROW/NCQL/  NELFM/NETA/  NLELEM/hiLETA/NOELEM/  NUETA/ JH(60)/ 

• )  KINBAS(122)/LA(122)/LE(502)/I l( 1000)/ IE (2000 ) 

CQMM0N/Gc.STLT/PPICE(130)/ICURX(130)/ISUMC,ITSINV/LISTL,  NBRAIC/ 

1  NS/NP/JFCQL(11),J»ROM(11),JP£L&M(11),MAXC(10),MAXC2(10) 

C 

IFKAS  =  0 
LI STL  =  0 
ISUMC  =  0 
INCVAL  =  INITBb 
C 

C  TEST  FOR  FATHOMING 
C 

ENTRY  8BENTR 
100  CALL  TESTX 

IF  (MSTAT  . F(i.  QBL)  GO  TO  200 
IF  (MSTAT  .EQ.  OL)  RETURN 

C  CURRENT  NODE  FATHOMED; BACKTRACK  TO  LAST  PROMISING  NODE  ON  LIST 
1*50  CALL  BKTRAK 

C  IF  LIST  IS  EMP1Y/  RETURN  TO  MAIN  (COMPUTATIONS  COMPLETED) 

IF  (LISTL  .EQ.  0)  RETURN 

C  USE  PRIMAL  SIMPLEX  METHOO  FOR  REOPTIMIZATION  AT  NEW  NODE 
CALL  NOOMAL(ITSINV) 

IF  (ITCNT  .CE.  ITRFRQ)  GO  TO  2000 
GO  TO  100 

C  CURRENT  NODE  NOT  FATHOMED;  COMPUTE  PENALTIES 
C  BRANCHING  AT  CURRENT  NODE  IS  DONE  FROM  SUBROUTINE  PENLTS 
200  CALL  PENLTS 

IF  (IDIK)  400/150/400 

C 

C  kFINVFRT  CURkF.MT  BASIS 
1000  CALL  INVERT 
ITSTNV  =  0 
C 

C  DUAL  SIMPLEX  CYCLE 
C 

C  CHOOSE  PIVOT  ROW  IKOWP 
300  CALL  DCHUZR 

IF  (IROWP  .CT.  0)  GO  TO  400 
HSTAT  =  QPL 
GO  TO  100 

C  C  irnSE  PIVJT  COLUMN  JCOLP 
-40C  CALL  DCriUZC 

IF  (JCOLP  .PU.  0)  GO  TO  150 

C  UPDATE  RIGHT-HAND  SIDES  TO  REFLECT  NEW  BASIS  RESULTING  FROM 
C  CUPRPNT  SIMPLEX  PIVOT 
CALL  HP PET A 
ITCET  =  ITCNT  ♦  1 
IF  (ITCNT  • GE.  ITRFRO)  GO  TO  2000 
ITSIHV  =  ITSINV  ♦  1 

Is  (NELT. M.GT.(  NENAX-NRQW)  .OR.  (  ITSTNV .GE. INVFRQ ) )  GO  TO  1000 
C  WRITE  JUT  NEW  bT A- VECTOR  FOP  CUPPFNT  SIMPLEX  PIVOT 
CALL  WHET A 
G1  TO  300 


i 
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C 

2000  LISTL  =  0 
RETURN 
F.ND 

- - 

SUBROUTINE  DCHUZR 

C 

C  SELECT  PIVOT  K3W  IROWP  F OR  CURRENT  DUAL- SIMPLEX  ITERATION. 

C  3FT  IROWP=0  IF  CURRENT  BASIS  IS  OPTIMAL.  OTHERWISE/  CHOOSE 

C  IROWP  TO  BE  THE  HOW  WITH  GREATEST  PRIMAL  INFEASIBILITY. 

C  SUBROUTINE  ADAPTED  PROF  INTEGER  PROGRAMMING  CODE  BB/  WRITTEN 

C  BY  GARY  A.  KOCHMAN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C 

IMPLICIT  REALM  (A/C/E-H/O/P/P-W/Z),  REALM  (B/D/X/Y)/ 

1  INTEGFF  M  (I-N/U) 

INTEGER  JM/KIVBAS/LA/LE/IA/IE 
I  NTF.G5K  IPAHI/INCUMB/IVBNO/IVID/IOBND 
DOUBLE  PRECISION  E(2000) 

RML  A(1000) 

C 

COMMUN/CCNSTS/  ZTOLZE/ZTOLPV/ZTCOST/ZTOLSM/NEGINF/ NEMAX/ NRNAX/QBL/ 
l  NLES/MTMAX/QA/QI/QF/QN/QSUB/QB/3C/Q£/Qti/QL/QQ/QR/QM,UG 
COMMON/ BBLIST/  DFPART(60) /RE VflND/ INCVAL/ I COL/ IV AL/ IDIR/ IP ART (122)/ 
1  INCUMP(130)/IYbJID(500)/IVID(500),TQ3ND(500)/NPIYOT/IPTYPE/IFEAS 
COMMON  XLB(122),XUB(122)/DE/DP/B(60)/X(60)/Y(60),YTEMP(60)/A/ 

1  E/MSTAT/IOEJ/IROWP/ITCNT/INYFRO/ITRFRQ/JCOLP/ 

2  NROW/NCOL/ NELEM/NETA/ NLELEM/ NLETA/ NUELEM/ NUETA/ JH(60)/ 

3  KINBAS(122)/LA(122)/LE(S02)/IA(1000)/IE(2000) 
C0MM0N/TIMEPS/1T0T/TSTURE/TIMELP/TIMEDR/TIMEDC/TIMINV 

C 

D  I  TIME  =  IHPTIM<1) 

IPOWP  =  0 
DP  =  -1.E10 
D)  1000  I=1/NR0W 

IF  (T  .EQ.  I OBJ)  GO  TO  1000 
ICOL  =  JH(I) 

IF  (X(I)  .LT.  (XLB( ICOL)  -  ZTCLZE))  GO  TO  100 
IF  (X(I)  .GT.  (XUB(ICOL)  ♦  ZTOLZE))  GO  TO  200 
GO  TO  1000 
C 

C  I  ASIC  VARIABLE  ON  ROW  I  FALLS  BELOW  ITS  LOWER  BOUND 

100  DE  =  XLB(ICOL)  -  X(I) 

IF  (D£  .LE.  DP)  GO  TO  1000 
IPTYPfc  =  0 
GO  TO  250 
C 

C  1 ASIC  VARIABLE  ON  ROW  I  EXCEEDS  ITS  UPPER  BOUND 

200  DE  s  XU)  -  XUB(ICOL) 

IF  (DE  .LE.  DP)  GO  TO  1000 
IPTYPE  =  -I 
C 

250  IROWP  =  I 

DP  *  DE 

1000  CONTINUE 

D  IT I ME 2  *  IHPTIM(I) 

0  TIMEDP  *  TIMOR  ♦  (  f  TIHE2- ITIME )/l 00000. 

RETURN 

END 


C 


SUBROUTINE  DCHU7.C 


nao  o  n  wiu  n  o  o  r:  n  n  r  n  n  n  n  n  n 
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SELECT  PIVOT  COLUMN  JCOLP  FOR  CURRENT  DUAL-SIMPLEX  ITERATION. 
SFT  JCOLP  -  0  IP  LP-PROBLFM  AT  CURRENT  NODE  IS  INFEASIBLE# 
OTHERWISE  CHOOSE  JCOLP  TO  MAINTAIN  PRIMAL-OPIIMALITT. 
SUBROUTINE  ADAPTED  FROM  INTEGER  PROGRAMMING  CODE  BB#  NRITTEN 
BY  GARY  A.  KOCHMAN  (OPERATIONS  RESEARCH#  STANFORD  UNIVERSITY) 

IMPLICIT  REALM  < A#C#E-H#0#P#R-W#Z)#  REALM  (8#0,X#Y)# 

1  INTPGERM  (I-N#Q) 

INTEGER  JH/NINBAS/LA#LE# IA# IE 
I  '>TEGER  IPART#  INCUMB#  I VB.ND#  I V  ID#  IOBND 
DOUBLE  PRECISION  F(2000> 

RFAL  A(1000) 

COMMON/CONSTS/  ZTOLZ£#ZTOLPV#ZTCOST#ZTOLSM#NEGINF#NEMAX#NRMAX#QBL# 
1  NL£S/NTMAX#QA#QI#QF#QN#QSUB#QB#QC#QE#QH#QL#QO#QR#QM#QG 
COMMON/BBLIST/  DFP ART(60)#  KEVBND# I NCV AL# ICOL# IVAL#IDIR# IP ART (122)# 
1  INCUHBI 130 )#IVBND(500 )# IVID(500># IOBND (500 )#NPI VOT# IPTYPE# IFEAS 
COMMON  XLB(122)#XTJB(122)#DF# DP#B(60)#X(60)#Y(60) #YTEMP(60 )# A# 

1  E#.*‘STAT#  109 J#  IROWP#  ITCNT#INVFRQ#  ITRFRQ#  JCOLP# 

2  NROW/NCOL#  NELEM/NETA# NLELEM#NLETA#NUELEM#NUETA#JH(60)# 

3  K INB AS( 122)#LA( 122 )#LF(502)# I A(1 000)# IE( 2000) 
C1MMGN/TIMERS/1T0T#TST0RE#TIMF.LP#TIMEDR#TIMEDC#TIMINY 

I TIME  =  IHPTIM(l) 

JCOLP  =  0 

IF  (IPTYPE  „EQ.  -1)  GO  TO  1000 

LEAVING  VARIABLE  FALLS  BELOW  ITS  LONER  BOUND;  COMPUTE  MAX  RATIO  DP 

DP  =  -1.E10 
m  500  J=l#NCOL 

IF  (KINBAS(J)  .GT.  0)  GO  TO  500 

IE  ((XUd(J)  -  XLB(J))  .LE.  ZTCLZE)  GO  TO  500 

F  =  J 

CALL  UNPACK(F) 

CALL  FTRAN(l) 

II  (KINBAS(J)  .EQ.  -1)  GO  TO  200 
IF  (Y(IROWP)  ♦  ZTOLPV)  225#225#500 
OC  IF  (Y( IROWP)  -  ZTOLPV)  500#225#225 

2S  Dc,  =  Y(IQBJ)/Y(  IROWP) 

IF  (PE  .LE.  DP)  GO  TO  500 
JCOLP  =  J 
PP  =  DF 

00  CONTI  *  UK 

STOCK  PIVOT  COL  JCOLP  IN  Y;  STORE  CHANGE  IN  INCOMING  VAR.  ICOL  IN  DP 
I*  (JCOLP  .EG.  0)  RETURN 
C ALL  UNPACK(JCOLP) 

CALL  FTRAN(l) 

1 0 JL  =  JH( IROWP) 

DP  =  (X(IRLWP)  -  XLB( ICOL) )/Y( IROWP) 

G 1  TO  2000 

LEAVING  VARIABLE  EXCEEDS  ITS  UPPER  ROUND;  COMPUTE  MIN  RATIO  DP 
1000  DP  *  1.E10 

r.n  1500  jai#NCnt 

IF  (FINBAS(J)  .GT.  0)  GO  TO  1500 

Ir  ( ( XU K ( J )  -  XLB( J) )  .LE.  ZTOLZE)  GO  TO  1500 


i 


K  =  J 

CALL  UNPACK(K) 

CALL  FTKAN(l) 

IF  (KINBAS(J)  .EQ.  -1)  GO  TC  1200 
T y  (Y(IKOWP)  -  ZTJLPV)  1500/1225/1225 
1200  IF  <Y(IfcOWP)  ♦  ZTOLPV)  1225/1225/1500 

1225  D t  =  Y(IOtJ)/Y(lRUWP) 

IF  (PE  .GE.  DP)  GO  TO  1000 

JCOLP  =  J 
DP  =  DE 

1 5CC  CONTINUE 

C 

r  STOKE  PIVOT  COL  JCOLP  IN  Y;  STORE  CHANGE  IN  INCOMING  VAR.  ICOL  IN  DP 
IF  (JCOLP  .EQ.  0)  RETURN 
CALL  UNPACK (JCOLP) 

CALL  FTRAN(l) 

ICOL  =  JH(IRGWP) 

UP  =  (X(IROWP)  -  XUB( ICOL) )/V ( I  ROMP) 

c 

2000  IF  (KIND AS (JCOLP)  .EQ.  0)  DE  =  DP  ♦  XLB(JCQLP) 

IF  (K I NBAS (JCOLP)  .EQ.  -1)  DE  =  DP  ♦  XUB(JCOLP) 

NPIVC1T  =  1 

D  IT  I  ME  2  =  IHPTIM(l) 

D  T IMF DC  =  TIMEDC  ♦  ( ITIME2-ITIME )/l 00000. 

PETUPN 

END 

C - 

SUBROUTINE  TESTX 

C 

C  TEST  LP-OPTIMAL  SOLUTION  AT  CURRENT  NODE  FOR  FATHOMING. 

C  FATHOMING  OCCURS  IF: 

C  (1)  LP  PROBLEM  AT  CURRENT  NODE  IS  INFEASIBLE  (MSTAT  =  QN); 

C  (2)  LP-OtT  OBJ.  VALUE  ♦  OBJ.  VALUE  FOR  PREVIOUS  PERIODS  ♦ 

C  OBJ.  ROUND  ON  SUCCEEDING  PERIODS  <=  OBJ.  OF  INCUMBENT 

C  (3)  LP-CPT  SUL.  SATISFIES  INTEGER  RESTRICTIONS  AND  NS  =  NP 

C  IF  THE  LP-CPT.  SOL.  IS  INTEGER  BUT  FATHOMING  DOES  NOT  OCCUR/ 

C  BRANCH  ON  (FIX)  ALL  NON  SLACK  VARS  AND  STORE  SUBPROBLEM  NS  IN 

C  PREPARATION  FOR  A  FORWARD  STEP.  SET  MSTAT  =  QE  TO  FLAG  THIS, 

f  SUBROUTINE  ADAPTED  FROM  INTEGER  PROGRAMMING  CODE  BB/  WRITTEN 

C  BY  GARY  A.  KOCHMAN  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 

C 

IMPLICIT  REALM  (A/C/E-H/O/P/R-W/Z)/  REALM  (B/D/X/Y)/ 

1  INTEGSRM  (I-N/Q) 

INTEGER  JH/KINBAS/LA/LE/IA/I5 
INTEGER  IP ART/ INCUMB/ IVBND/ IVID/ IOBND 
DOUBLE  PRECISION  E(2000) 

REAL  A(1000) 

C 

COMMUN/CGNSTS/  ZTOLZE/ZTOLPV/ ZTCQST/ZTOLSM/NEGINF /NEM AX / NRMAX/QBL/ 
1  HLKS/NTMAX/QA/QI/QF/UN/QSUB/QB/QC/QE/QH/QL/QO/QR/QM/QG 
COMHON/BHLIST/  DEPART  (60  >,-  RKVBND,  INCVAL/ ICOL/  IV  AL/ IDIR/  IP  ARTC  122)/ 
1  INCUMB(130)/IVBND(500)/IVID(500)/IONND(500)/NPIVOT/IPTTPE/IFEAS 
COMMON  XLB(12?)/XUB(122)/D£/DP/B(60)/X(60)/Y(60)/YTENP(60)/As 

1  E/ MS TAT/ IOBJ/ IROWP/ ITCNT/ INVFR Q/ IT RFRQ/ JCOLP/ 

2  NROW/NCOL/NELFM/NFT A/NLELEM/NLET A/NUELEM/NUETA/ JH(60)/ 

3  KINBAS(122)/ LA( 122)/LE( 502)/ I A( 1000)/ IE (2000) 

COMMON/ GESTLT/ PRICE ( 130)/ ICURX( 130)/ ISUMC/ITSINV/LISTL/ NBR ARC/ 

1  MS/NP/ JFCOL( 1 1 )/ JFPOW( 11 )/ JFELEM( 11 )/MAXC(10)/ MAXC2( 10) 
COMMON/TIMERS/ ITQT/TSTORE/ TTMELP/TINEDR/TIMEDC/TINIRV 
C 


non  r.  r.  *.  »i  r.  t~.  «-»  •j'  n  n  <"»  c 
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C  INCUMB  =  BEST  SOLUTION  POUND  SO  FAR;  INCVAL  IS  ITS  OBJECTIVE  VALUE 

C  ICURX  CONTAINS  CURRENT  VALUES  FOR  VARIABLES  IN  SUBS.  1,...,NS-1 

C  OFPART(I)  L  IFART(  1)3  IS  THF,  FRACTIONAL  C  INTEGRALS  PART  OF  X(I) 

C  1VAL  =  LP-UPT.  OBJECTIVE  VALUE  FUR  SUBPROBLEM  NS  ON  CURRENT  BRANCH 

C  J90UND  =  BOUND  ON  MAX.  OBJECTIVE  VALUE  (USING  MAXC) 

C  1SUMC  =  OBJECTIVE  VALUE  FOR  SUBPROBLEMS  1,...,NS-1  ON  THIS  BRANCH 
C  I^F  A3  =  1  IFF  AN  INTEGRAL  SOL.  Ml  Til  INCVAL  >  INITBD  HAS  BEEN  POUND 

C 

res T  FOR  FATHOMING  IN  MAPS  (1)  AND  (2) 

IF  (MSTAT  . L'Q.  ON)  GO  TO  2000 
DP  =  X(IOBJ)  ♦  ZTOLZE 
IVAL  =  101 N^ (DP ) 

IF  (DP  .LT.  0.)  1 VAL  *  IVAL  -  1 

IF  (( MAXC2 (NS) +IVAL+ ISUMC)  .LE.  INCVAL)  GO  TO  2000 
IF  (NS.EO.NP)  GO  TO  50 
JBQUND  =  IVAL  ♦  MAXC(NS) 

Ic  (JBiJUND.LT. 0)  JBOUND  =  JBOOND  -  1 
J  UJUND  =  J  BOUND  /  2 

IF  ((JBOUND  *  ISUMC)  .LE.  INCVAL)  GO  TO  2000 

COMMUTE  INTEGER  AND  FRACTIONAL  PARTS  OF  EACH  BASIC  VAR. 

0  100  I=l,VROW 

IP ART ( I )  =  IOINT(X(I)  ♦  ZTOLZE) 

IF  (X(I)  .LT.  -ZTOLZE)  IPART(I)  =  IPART(I)  -  I 
MTEMP1  =  IP ART (I ) 

CC  DFPART(l)  *  X(l)  -  FLOAT (NTF.MP1) 

CHECK  FOP  ALL-INTEGER  SOLUTION 

DO  200  1*1, "ROM 

TF  ( JH ( 1 )  .LE.  NROM)  GO  TO  200 
IF  (DFPART(I)  .GE.  ZTOLZE)  RETURN 
0(y  CONTINUE 

SOLUTION  ALL-TNTEGcR:  CHECK  FOR  COMPLETE  SOLUTION 
Ip  (KS.LT.NP)  CO  TO  400 

*IKN  IMPROVED  INTEGER  SOLUTION  TO  ALL  PERIODS  REACHED. 

OUTPUT  OPJ.  VAL.  AND  COMPUTATION  TIME  REQUIRED  TO  REACH  IT. 

0 

INCVAL  =  IVAL  ♦  ISUMC 
JTIME  =  I'lP^IMd  ) 

TOPT  =  (JTIML-IT0T)/100000. 

IF  ( IPE AS. F.Q.O)  «RITE(21/  1 ) 

1  FORMAT  (*  INTERMEDIATE  SOLUTIONS  POUND') 

WRIT"  (21/2)  TOPT,NBRANC, INCVAL 

K1KMAT  (*  TIME  =',F7.2, '  SECONDS;  BRANCHFS  *',110,';  INCVAL  *', 
l  110) 

IFEA3  =  1 

C  STQnr  NSW  INCUMBENT  SOLUTION 
K  =  JFCOL(NS)  -  1 
00  300  J*l,V 

JOU  INCUMii(J)  =  ICURX(J) 

in  350  J*l,"COL 

IF  (riNBAS(J))  320,330,353 
120  INCUMB(K4J)  a  TD1NT( XUB( J) ) 

GJ  TO  350 
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330  TNCUMl,(K*J)  =  ID  INK  XLB(.J)  ) 

350  CONTINUE 

00  3o0  T=1,NK0W 

ICOL  =  JH( I ) 

360  r»CUMJ(K*ICOL)  =  IPAPT(l) 

Cn  TO  2000 

C 

C  \  fMKTIAL  INTEGFk  SOLUTION  HAS  BEEN  (CACHED.  BRANCH  ON  ALL 
c  JO'/SLACK,  (ISFIXr.n  VARIABLES  AND  SAVE  SOLUTION  IN  ICURX. 

c 

400  *  =  JF:nL(WS)  -  1 

r  tiRAVCK  ON  ILL  f.ONRASIC  NOVSLACK  VARIABLES  FIRST 
DO  500  J=l,NCOL 

Ir  (KINBAS(J))  420,440,500 
r  VAPlAJLi-  NIN’USIC  AT  UPPER  BOUND;  BRANCH  UP 

42C  ICUVX(IUJ)  =  IDINT(XUB(J)> 

It  (J.LE.NROW  .OR.  (XUR(J)-XLB(J)). LE.ZTOLZE)  GO  TO  500 
ID  IK  =  -1 

R-.VbND  =  SNGL(  XUB(  J)  ) 
r.)  TO  4b0 

C  VARIABLE  NONB AS IC  AT  LOWER  BOUND;  BRANCH  DORN 

44C  ICURX(K*J)  =  IDI NT(XLB( J) ) 

IF  (J.LE.NROW  .OR.  (XUB(J)-XLB(J)) .LE.ZTOLZE)  GO  TO  500 
I DIR  -  1 

RSVBNO  =  SNGL( XLB( J) ) 

460  ICOL  =  J 

CALL  BRANCH 
500  CONTINUE 

C 

C  STURE  AND  BRANCH  ON  ALL  NONSLACK  BASIC  VARIABLES 
DO  600  1=1, NROW 

ICOL  =  JH( I ) 

ICURX(X*ICOL)  =  IPAKT(I) 

IF  ( TCOL.LE.NKOi)  GO  TO  600 

IF  ((XUB(ICOL)-XLB(ICOL)). LE.ZTOLZE)  GO  TQ  600 
IF  (<X(I)-XL9(ICOL)). LE.ZTOLZE)  GO  TO  520 
IF  ((XUB(ICQL)-X(I) ) .GT. ZTOLZE)  GO  TO  550 
C  VARIABLE  BASIC  AT  UPPER  BOUND;  BRANCH  UP 

IDIfi  =  -1 

RSVBND  =  SNGL(XUB(ICOL)) 

GU  TO  590 

C  VARIABLE  BASTC  AT  LONER  BOUND;  BRANCH  DOWN 

52C  IDIR  =  1 

REVB?!D  =  SN*>L(  XLR(  I  COL  )  ) 

GO  TO  590 

C  VARIABLE  BASTC  BETWEEN  BOUNDS;  BRANCH  TO  FIX  IT 

550  IDIR  =  -1 

RKVBND  =  FLOAT(IPARTd)) 

CALL  BRANCH 
IDIR  =  1 

590  CALL  BRANCH 

600  CONTINUE 

C  STORE  CURRENT  SUBPROBLEN  AND  THE  CURRENT  OBJ.  VALUE 
CALL  STORE 
IS1IMC  =  ISUPC  ♦  I VAL 
MSTAT  =  QE 
kOTUKN 
C 

C  CURRENT  PRQBLFfc  f»n  LUNGER  OF  INTEREST 
C 


n  r.  o  c.  ci  -■>  o  t~,  <~i  n  o 


2000  MSTAT  =  01 
RETURN 
E»'0 

SUBROUTINE  PENLTS 

COMPUTE  TOMLIN'S  IMPROVED  OP*  AND  DOWN-  PENALTIES  AND  THE 
GOMORY  PENALTY  FOR  EACH  NONINTEGER  BASIC  VARIABLE.  THEN  CHECK 
FOR  FORCED  BRANCHES  ON  BOTH  BASIC  AND  NONBASIC  VARIABLES.  IN 
THE  ABSENCE  OF  FORCED  BRANCHES  ON  BASIC  VARIABLES/  ADD  TO 
EACH  PENALTY  THE  PRICE  OF  THE  CORRESPONDING  OFFDI AGONAL  COLUMN 
OF  THE  NEXT  SUBPROBLEM  (NS  ♦  1).  THEN  CHOOSE  AS  BRANCHING 
VARIABLE  THE  ONE  WITH  LARGEST  ASSOCIATED  UP-  OR  DOWN-PENALTY. 
TAKE  THE  FORWARD  BRANCH  IN  THF  DIRECTION  OPPOSITE  TO  THIS 
MAXIMUM  PENALTY  WHILE  ADDING  THAT  VARIABLE  TO  THE  LIST  WITH 
THE  APPROPRIATE  BRANCH  DIRECTION  (IVID)  AND  BOUND  (IOBND). 

THE  BRANCHING  PROCESS  ITSELF  IS  DONE  IN  SUBROUTINE  BRANCH. 
SUBROUTINE  ADAPTED  FROM  INTEGER  PROGRAMMING  CODE  BB/  WRITTEN 
C  BY  GARY  A.  KGCHMAN  (OPERATIONS  RESEARCH#  STANFORD  UNIVERSITY) 

l 

IMPLICIT  REALM  ( A/C/E-H/O/P/R-W/Z)/  REALM  (B/D,X,Y)/ 

1  INTEG.ERM  (I-N/Q) 

INTECEK  JH/KINBAS/LA/LE/IA/IE 
I NTEGER  IPAPT/ INCUMfi/ IVBND/IVID/ IOBND 
DOUBLE  PKECTSIGN  F(2000) 

HEAL  A(1000) 

REAL  DU(60)/PO(60)/PG(60) 

LOGICAL  FORCED 

C 

COMMON/CONSTS/  ZT0LZE/ZT0LPV/ZTC0ST/ZTULSM/NEGINF/NEHAX/NRHAX/Q8L/ 

1  i-LES/NTMAX/OA/GI/QF/QN/QSUB/QB/QC/QE/QH/QL/QO/QR/QM/QG 
C  Jl  MON/PBLIST/  DFPART(60)/P.EVBND/INCVAL/  ICOL/IV  AL/1DIR/  IP  ART  (122  )/ 

1  INCUML(13C)/IVEND(500)/ IVID(500)/ IOBND<500)/NPIVQT/IPTYPE/IFEAS 
COMMON  XLB(122)/XUB(122)/0E/DP/B(60)/X(60)/Y(60)/YTEMP(60)/A/ 

1  h/MSTAT/ IOLJ/ TROWP/ ITCNT/ INYFPa/ ITRFRQ/ JCOLP/ 

2  JROW/NCOI./NF.LFM/NFTA/NLELFM/NLET  A/NDELEM/NUETA/ JH(60)/ 

0  K.I'JoAS(122)/LA(122)/LE(S02)/IA(1000)  /  IE  (2000  ) 

Cf)!iMUN/CESTLT/FRICE(130)/ICURX(  130)/ ISUMC/ITS INV/LISTL/  NBRANC/ 
i  RS/NP/JFCOL(11)/JFROW(11)/JFELEM(1 1)/MAXC(10)/MAXC2(10) 
c 

C  <’J(T)/PD(I)/Pr(I)  ARE  THE  UP/DUWN/  AND  GUMORY  PENALTIES  FOR  VAR  JH(I) 

C  .-'OVDHD  =  REVISED  BOUND  ON  THE  BRANCH  VAP  ICOL?  IT  IS  PASSED  TO  BRANCH 

:  l'  I0IR  =  -1  BRANCH  (JO/  *1  BRANCH  DOWN/  0  FATHOMING  HAS  OCCURRED 

C  1/AL  =  LP-OPT.  OtiJfcCriVE  VALUE  FUR  SUBPOU8LEM  NS  (JN  CURRENT  BRANCH 
C  JPOUnD  -  BOUND  ON  MAX.  OBJECTIVE  VALUE  FOR  SUBPROBLEMS  NS/.../NP 
r  ISU'C  =  JBJECT1VE  VALUE  FOR  3UJPK0*LEMS  1/.../NS-1  ON  THIS  BRANCH 
C  /HCTJk  Y  CONTAINS  APPRUPRIATE  COLUMN  OF  CURRENT  SIMPLEX  TABLEAU 
C 

UD  10  I =1/ N ROW 

IF  (DFPART(I)  .LT.  ZTDLZE)  CO  ?U  5 
PU(I)  s  1.E6 
PU(I)  *  1.F6 
PG( I )  =  1.E6 

go  in  10 

G  PU(I)  =  o. 

PD( I )  =  C. 

PG( I )  s  0. 
in  CONTINUE 

r 

r  *11”  LOOP:  CALCULATE  PENALTIES  FOP  EACH  ELICIRLE  INCOMING  VAN.  J 
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C 

i'J  1000  J=1,NCUL 

II  (MNBAS(J)  .GT.  0)  GO  TO  1000 

JF  ((XUb(J)  -  XL’’(J) )  .LE.  ZTOLZE)  GO  TO  1000 

^  -  J 

CALL  UNPACK(K) 

CALL  FT.UN(l) 

IF  (KINBAS(J)  •  EQ.  0)  GQ  TO  30 
O')  20  I=l,NRO* 

?0  V ( I )  =  -Y(I) 

C 

C  C.IF.CK  FOR  FORCED  BRANCH  UN  NONBASIC  VARIABLE  J 

C 

JO  IF  (J  .LE.  NROV)  GO  TO  60 

DP  =  X(iOBJ)  -  Y(IObJ)  ♦  ZTOLZE 
IVAL  =  lDIMT(UP) 

Ir  (DF  .LT.  0.)  IVAL  =  IVAL  -  1 
IF  ((MXC2(NS)  +  1VAL  +  ISUHC)  .LE.  INCVAL)  GQ  TO  50 
IF  (NS.EO.NP)  GO  TO  60 
JBOUNP  =  IVAL  ♦  MAXC(NS) 

IF  (JfcOUND.LT.O)  JROUND  =  JBOUNO  -  1 
JBOUNO  =  JBOUND  /  2 

IK  ((JBOUND  ♦  ISUMC)  .GT.  INCVAL)  GO  TO  60 
r-0  IDIK  =  2*KINRAS(J)  ♦  1 

IF  (IUIR  .EU.  -1)  REVBNO  =  SNGL( XUB( J) ) 

It  (TUIK  .EQ.  1)  REVBND  =  SNGL(XLB(J)) 

I COL  =  J 
CALL  BRANCH 
GC  TO  1000 


C  COMPUTE  PENALTIES  ON  BASIC  VARIABLES  JH(I>,  FOR  ICOL  s  J 

*>0  DO  500  I=1,NRQW 

IF  (JK(I)  .Lh.  NROW )  GO  TO  500 
IF  (DFPAKT(I)  .LT.  7.TQL7E)  GO  TO  500 
C  COMPUTE  UP  PENALTY  FOR  JH(T) 

100  IF  ( Y( I )  .GT.  -ZTOLPV)  GO  TO  200 

DE  =  Y(IOBJ)*(DFPART( I)  -  l.)/Y(I) 

IF  (DE  .LT.  Y ( IOBJ ) )  DE  =  Y(IOBJ) 

IF  (DF  .LT.  PU(I))  PU(I)  =  DE 
GO  TO  300 

COMPUTE  DOWN  PENALTY  FOR  JH(I) 

200  IF  (Y(I)  .LT.  ZTULPV)  GO  TO  300 

DE  =  Y ( IOBJ )*OFP APT ( I )/Y(T  ) 

IF  (DF  .LT.  Y(IOBJ))  DE  =  Y(IOBJ) 

IF  (DE  .LT.  PD(I))  PD( I )  =  DE 
r  COMPUTE  GOMORY  PENALTY  FOR  JH(I) 

300  DP  =  DABS (Y(I ) ) 

IF  (DP  .LE.  ZTOLZE)  GO  TC  500 
NTEMP1  =  IDINT(DP) 

DP  =  DP  -  FLUAT(NTEMPl) 

IF  ((DP  .GT.  ZTOLZE)  .AND.  (DP  .LT.  l.-ZTOLZE))  GO  TO  330 
IF  (J.NE.I)  GO  TO  500 
IF  (Y(I)  .LT.  0.)  GO  TO  320 
OK  =  Y (  IllBJ) *DFPAPT( I )  /  Y(I) 


CO  TO  350 

370 

DE 

=  Y ( I OBJ )* ( 1 .  -  DFPART(I))  / 

GO  TO  350 

330 

IF 

(YU)  .LT.  0.)  DP  s  l.  -  DP 

IF  (DP  .GT.  DFPART(I))  CO  TO  340 
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UE  =  Y(TOBJ)*DFPART(I)/DP 


GO  TO  350 

340 

DE  s  Y(  I(JBJ)*(  1 .  - 

DFPART(I))/(1.  -  DP) 

350 

IF  (DE  .LT.  PC(I)) 

PG( I )  =  DE 

500 

CONTINUE 

10CC 

CONTINUE 

C  COMPUTE  LARGEST  CONOR Y  PENALTY  AND  TEST  FOR  PATHOMING 
C 

PEN  =  0. 

no  2000  I=l,NROW 

IF  (JH(l)  .LE.  NROM)  GO  TO  2000 
IF  (PG(I)  .GT.  PEN)  PEN  =  PGQ) 

1000  CONTINUE 

OP  =  X(IOBJ)  -  PEN  ♦  ZTOLZE 
IVAL  =  IL)INT(DP) 

IF  (OP  .LT.  0.)  IVAL  =  IVAL  -  1 

1?  ((MAXC2(NS)*I VAL+ISUMC)  .LE.  INCVAL)  GO  TO  2050 
If  (N3.Eu.NP)  GO  TO  3000 
JliOUND  =  IVAL  ♦  MAXC(NS) 

If  (JROUNO.LT.O)  JfcOUND  =  JBOUND  -  1 
JIO'JND  =  JEDUMi  /  2 

ie  ( ( jbuunc  ♦  isumc)  .gt.  incval)  go  to  3000 

^05 0  iniR  *  0 
RETURN 

C  PROBLEM  NOT  FATHOMED:  CHECK  FOR  FORCED  BRANCHES  ON  BASIC  X(I) 
r 

1000  FORCED  =  .FALSt. 

DO  3*300  I=1,NRUN 

IF  (JH(I).LE.NROW  .OR.  DFPART( I) . LF. ZTOLZE)  GO  TO  3900 
IF  (PU(I)  . GT.  PD( I ))  GO  TO  3600 
0°  =  X(10flj)  -  PD( I )  ♦  ZTOLZE 
MTEM.P1  *  IDINT(DP) 

IF  (Dr  .LT.  0.)  NTEMP1  =  NTe«Pl  -  1 

IF  ((MXC2(MS)«-NTEMP1*ISUMC)  .LE.  INCVAL)  GO  TO  3050 
IF  (MS.EQ.NP)  GO  TO  390C 
JROUHD  =  NTEMP1  ♦  BAXC<NS) 

IF  (JhOUNU.LT.O)  jeoUND  =  JBOUND  -  1 
JHOON'D  =  JBUUND  /  2 

IF  ( (JBUUND  ♦  ISUMC)  .GT.  ThiCVAL)  GO  TO  3900 
f  1RCED  3RANCH  UP  ON  r<I) 
jO'jO  TVAL  =  NTcMPl 

II)  II:  =  *1 

N1FMP1  =  IPART(l)  ♦  1 
G  i  TO  3700 
C 

360*  DP  =  A(IOBJ)  -  PU( T)  ♦  ZTOLZE 

NTEMU  *  IDINTCDP) 

I!-  (Df  .LT.  0.)  NTEMPl  =  NTEMFl  -  1 

IF  ((►AXC2(NS)»NTLMPl>tSUM,C)  .LE.  INCVAL)  GO  TO  3650 
Tr  (»'5.hG.MP)  GO  TO  3900 
JrfflUVD  *  {JTEMPl  ♦  MAXC(NS) 

IF  ( JFOUVD.LT.0)  JDJUND  =  JbOUND  -  1 
.1  OUMt,  =  JBOUND  /  2 

TF  ((J°nWND  ♦  ISUMC)  .GT.  INCVAL)  GO  TO  3900 
C  FtIRCF'J  dRHT'i  DOWN  ON  X(I) 

J5'r-  IVAL  a  NTfeHi-i 

Ti'.IR  *  1 

NTFMPI  a  IPAPTU) 
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3700  IFOWP  =  T 

l  COL  =  JH(IROWP) 

Ri  VBV'O  =  FLOAT(NTEMPl) 

FORCED  =  .TRUE. 

CALI.  BRANCH 
3900  CONTINUE 

IF  (FORCED)  GO  TG  5000 
C 

C  MO  FORCED  BRANCHES!  CHOOSE  BRANCHING  VAR.  AND  DIRECTION 
C 

PPN  =  0. 

nuWP  =  0 

C  DETERMINE  BASIC  VARIABLE  JHdROWP)  WITH  MAX.  UP-  OR  DOWN-PENALTY, 

C  ADDING  IN  A  PENALTY  FROM  THE  NEXT  SUBP30RLEM  (STORED  IN  PRICE) 

DO  4900  1=1, BROW 

IF  (JH(  I)  .LE.I1ROW  .OR.  DFPART( t) .LF.ZTOLZE)  GO  TO  4900 
UPPEN  =  PU(I)  ♦  PRICE(JFCOL(NS)  ♦  JH(I)  -  1) 

IF  (UPPEN  .CT.  PD(I))  GO  TO  4600 
IF  (PO( I )  .LE.  PEN)  GO  TO  4900 
PEN  =  PD( I ) 

IKOWP  =  I 
1 0 1 R  =  -1 

N TEMPI  =  IPART(I)  ♦  1 
REVBND  =  FLOAT(NTEMPl) 

GO  TO  4900 

46C0  IF  (OPPEN  .LF.  PEN)  GQ  TO  4900 

PEN  =  UPPEN 
TkOWP  =  I 
TDIR  =  1 

NIEMP1  =  IPAPT(I) 

Rt  VB^'O  =  FLOAT(NTEMPl) 

49C0  CONTINUE 

IF  (IRGWP  .CT.  0)  GO  TO  4950 

C  ,.AC'«  UP-  AND  OUWN-PENALTY  =  0.  (DUAL-DEGENERACY)  CHOOSE  ANY 
C  NONINTEGER  BASIC  VARIABLE  AS  BRANCHING  VARIABLE  ICOL 
DO  4910  IRQWP=1, NRUW 

IF  (JH(IRONP)  .LE.  NROW)  GO  TO  4910 
IF  (DFPART(IRO.IP)  .GE.  ZTOLZE)  GO  TO  4920 
4910  CONTINUE 

4920  ID IR  =  1 

N TEMPI  =  TPART(TROWP) 
h  KVB'lD  =  FLUAT(NTEMPl) 

4950  IF  ( ID IR  • EQ.  1)  PEN  =  PU(IRONP) 

ICllL  =  JH(IPOWP) 

DP  =  X(IOBJ)  -  PEN  ♦  ZTOLZE 
NTEMP1  =  IDINT(DP) 

IF  (DP  .LT.  0.)  NTEMP1  =  NTEMP1  -  1 
IF  ( I VAL  .GT.  NTEMP1)  IVAL  =  NTEMP1 
C  BRANCH  ON  CHOSEN  VARIABLE 
CALL  BRANCH 

5000  IF  (IDIR  .LO.  -1)  IPTVPE  =  0 
IF  (IDIR  . EQ.  1)  IPTYPE  =  -1 
RETURN 
END 

SUBROUTINE  BRANCH 
L 

C  BRANCH  ON  VARIA3LF  ICOL  AS  DPTFRM I  NED  IN  SUBROUTINE  PEKLTS 

C  SUBROUTINE  ADAPTED  FROM  INTEGER  PROGRAMMING  CODE  BB,  BRITTEN 

r  BY  GARY  A.  rOCHMAN  (OPERATIONS  RESEARCH,  STANFORD  UNItERSITT) 


i  n  on  r .  n  o  n  o  n  n  o 
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IMPLICIT  RE*L*4  (A/C/E-H/O/P/R-K/7.  )/  »EAL*8  (B/D/X/Y)/ 
1  1MTEGERM  (I-N/Q) 

INTEGER  JH/ K INBAS/LA/LE/ I A/TE 
INTEGER  I  PART/ IN CUM 8, IVBND/IVID/IQBND 
DOUBLE  PRECTSIUN  £(2000) 

HEAL  A(1000) 


j  n  "’u  n  /  »ouio  a/  ur  rnn  i  \  uv/#ivci  vo  t  x  vub/  x  v  Mb/  ivin/ 

lNCUHB(130)/IVBND(500)/IVID(S00)/IOBND(S00)/NPIVOT/IP 
1MMQN  XL 8(122)/ XUh( 1 22)/ DE/ DP/ 8( 60)/ X(60 )/ Y(60 )/ YTEMP( 
E/ MS TAT/  10UJ/  TROWP/  I  TON T/  I  .)VFRQ/  ITRPRQ/JCQLP/ 
RROW/NCOL/NFLEH/ NET A/ NLFLEM/ MLETA/NUELEM/NUETA/ JH(60 ) 
KIN»>AS(122)/LA(122)/LE(502)/IA(1000)/IE(2000) 


ICOL  INDEXES  RRANCtiING  VARIABLE  CHUSEN 

IPIH  =  -1  MEAN'S  HRANCH  UP/  =  *1  MEANS  BRANCH  DOWN 

IVID  STORES  BRANCH  VARIABLE/  OPPOSITE  DIRECTION/  AND  SUBPROBLEN  # 

1'Jb'lD  STORES  AN  QHJECTl VE  FUNCTION  BOUND  ON  THE  OPPOSITE  BRANCH 

IV8SU  STOkKS  VARIABLE  BOUND  XUri  OR  XLB  FOR  OTHER  BRANCH  DIRECTION 

AJD  OPPOSITE  DIRECTION  TO  LIST 
NURANC  =  NBP  ANC  ♦  1 
L13TL  =  LI5TL  ♦  1 

IF  (IDIR  .EG.  -1)  IVBND(LISTL)  =  IDINT(XLB( ICOL)  ♦  ZTOLZE) 

IF  (IDIR  .EQ.  1)  TVPND(LISTL)  =  IDINT(XUB( ICOL)  ♦  ZTOLZE) 
miHLlSTL)  =  IDIR  *  (ICOL  «■  (NS  *  1000)) 

IfJbNO(LISTL)  =  I VAL 

r'SVISF  BOUNDS  ON  BRANCHING  VARIABLE  FOR  FORWARD  DIRECTION 
I*  (IPIP  .EQ.  -I)  XLB(ICUL)  =  DBLE(REVBND) 

Ip  (I11IH  .EQ.  1)  XUB(ICOL)  =  08LE(  REVBND ) 

KKTUV‘1 

END 


SUBROUTINE  RKTRAK 

BACKTRACK  TO  A  PROMISING  (IMFATHOMED)  NODE  FROM  THE  LIST  OF 
S^UVrcD  NODES.  EMPLOYS  LAST- IN-FIRST-UUT  (LIFO)  SELECTION  RULE 
SUF!HOUTI*’E  AD AnTSD  FROM  INTEGER  PROGRAMMING  CODE  BB/  WRITTEN 
UY  GARY  A,  KUCHMA*)  (OPERATIONS  RESEARCH/  STANFORD  UNIVERSITY) 


IMPLICIT  RfcALM  (A, 
I*!T FGEi\*4  (I-N/Q) 


,C/E-H/0/P/R-W/Z)/  RbAL*8  (9/D/X/1 
) 


1  I*!T FGEi\*4  (I-N/Q) 

1 N  TEGDr)  JH/fiNd*S/LA/LK/IA/IE 
I NTFGeV  IPAPT, INCUMB/IVBND/1 VID/IOB 
ipT-PLE  P  < EC  I  SION  S(2000) 
h'AL  A<  1000 ) 
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3  KIN!H<;(122),HU2?),LF.<502)/I»(1GOQ>/IE(2000> 

C)MMON/GESTLT/PRICii(130),ICURX(130),ISUMC,  ITSI NV/LISTL, NBR ANC, 
1  HS,NP,JFCOL(11),JFROW(11),JFELEM(11),MAXC(10),MAXC2<10> 

r 

N TEMP 3  =  0 

C  IF  LIST  IS  EMPIY,  RETURN  ( COMPUTATIONS  COMPLETED ) 

SO  I*  (LI STL  .KQ.  0)  RETURN 

C  LET  NEXT  NODE  FROM  LIST.  CHECK  ITS  SURPRQ8LEM  NUM3ER. 

TOOL  =  lVTD(LISTL) 
nunuo  =  iaes(icol)/iooo 

IF  (1SURN0.FG.RS)  GO  TO  7) 

C  .  E  \ RF.  BACKTRACKING  TO  X  PREVIOUS  SURPROBLEM 
T"»  =  ISUtJHO 
C ILL  RESTOR(O) 

70  IF  ((*AXC2(NS)*I031D(LISTL)*ISUMC)  .LE.  INCVAL)  GO  TO  2000 
IF  (.iS.E7.RP)  GO  TO  RO 
JHQUUi)  =  lOBND(LISTL)  ♦  MAXC(MS) 

IF  (JROUND.LT. 0)  JriOU N J  =  J3JUND  -  l 
JROUND  =  JBOUND  /  2 

IF  ( (JBOUND  ♦  1SUMC)  .LE.  INCVAL)  GO  TO  2000 
90  IF  (ICOL  .LT.  0)  GO  TO  100 
I COL  =  ICOL  -  (NS  *  1000) 

c 

C  liPAVCH  DIRECTION  WAS  DOWN.  RESTORE  UP  DIRECTION  BOUNDS. 

NTEKP1  =  IDINT(XLR(ICQL)  «•  ZTOLZE) 

NTEMP2  =  IVSKD(LISTL) 

XLP(ICOL)  =  XUL'(ICOL)  ♦  1. 

XUB(ICOL)  *  FL0AT(NTEMP2) 

IF  (KINRAS (ICOL)  .GT.  0)  GO  TO  1000 
KINPAS(ICOL)  =  0 
NTEMP3  =  1 
CO  TO  1000 
C 

r  BRANCH  DIRECTION  WAS  UP.  RESTORE  LOWER  DIRECTION  RJUNDS. 

IOC  ICOL  =  -  (ICOL  ♦  (NS  *  1000)) 

STfcMPl  =  IDINT(XUB(ICOL)  ♦  ZTOLZE) 

N TEMP 2  =  IVBNO(LISTL) 

X 7 B (ICOL)  =  XLB(ICOL)  -  1. 

XLa(ICOL)  =  FLOAT(NTEMP2) 

Ip  (KINPAS(TCOL)  .GT.  0)  GO  TO  1000 
KIN” AS( TOOL)  =  -1 
HTEMP3  =  1 
C 

C  MARK  OLD  BRANCH  AS  FATHOMED 

C 

1000  TVID(LISTL)  =  -IVID(LISTL) 

IVBVO(LISTL)  s  NTEMP1 
I OBND(LISTL)  =  NECINF 
C 

C  UPDATE  X  IF  NECESSARY 

C 

IF  (NTEMP3  .RE.  0)  CALL  UPDATX 
RETURN 
C 

C  NODF  FATHOMED:  UPDATE  VAR.  ROUNDS  AND  BACKTRACK  AGAIN 

C 

2000  IF  (ICOL  .LT.  0)  GO  TO  2100 
ICuL  *  ICOL  -  (NS  *  1000) 

N TEMPI  =  I VRhD(LISTL) 

IF  (KINBAS(ICOL))  2010,2050/2050 
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2010  *1  rtWP3  *  1 

OP  =  X'JB(ICOL)  -  XLB(ICOL) 

DY  =  FLOAT(NfEMPl)  -  XUB(ICOL) 

IF  (DP  .LT.  OY)  KINBAS(ICOL)  =  0 
2C50  X'IB(ICOL)  =  FLOAT(NTEMPl) 

GO  Til  3000 
C 

2100  ICOL  =  -(ICOL  ♦  (NS  *  1000)) 

NTF.MP1  =  I  VRNO(LISTL) 

IP  (KINB  AS(TCQL) )  2150/ 2110/  2150 
2110  NTFMP3  =  1 

D7  =  XLR(ICOL)  -  FLOAT(NTEMPl) 

DP  =  XUB( ICOL)  -  XL3(IC0L) 

IP  (DP  .LT.  DY)  KINBAS(ICOL)  *  -1 
2150  XLB(ICOL)  *  FLOAT( NTEMP1) 

C  CONTINUE  BACKTRACKING 
JOOO  LISTL  =  LISTL  -  1 
on  TO  150 
F.NL) 

SUBROUTINE  vpupup 
C 

r  OUTPUT  OPTIMAL  SOLUTION  AMD  CORRESPONDING  OBJECTIVE  VALUE. 

C  SUBROUTINE  ADAPTED  FROM  INTECER  PROGRAMMING  CODE  BB/  WRITTEN 

C  BY  GARY  A.  KOC11MAN  (OPERATIONS  RESEARCH/  STANPORD  UNIVERSITY) 

C 

IMPLICIT  KF.ALM  (A/C/E-H/Q/P/P-W/Z),  REAL*8  (3/0/X/Y)/ 

1  INTEGER**  (I-N/Q) 

INTEGER  jH/KINfcAS/LA/LE/IA/IE 
INTEGFR  IPART/INCUMB/IVBND/IVID/IOBND 
DOUBLE  PRECISION  F(2000) 

REAL  »(1000) 

COMMON/ BBL1 ST/  DFP ART(60) /RE VBND/ I NCV AL/ ICOL/ IY AL/ ID IR/ IP ART( 122)/ 
1  INCUMo(13C)/TVBND(500)/IVID(500)/T0BND(500)/NPIV0T/ IP TYPE/ I FE AS 

COMMON  XLO(122)/XUb(122)/Dfc;/OP/B(60)/X(60)/Y(60),YTEMP(60)/A/ 

1  E/MSTAT/ lOtJ/IKOWP/ITCNT/INVFPQ/ITRFRQ/JCOLP/ 

2  NROV/*1COL/NELFM/  NETA/NLELEN/NLET  A/NUELEM/NUETA/  JH(60)/ 

I  K INS  AS( 1 22)/LA (1 22)/LE(502)/ 1 A( 1000)/ IE ( 2000) 

C'JMMON/GESTLT/PRlCE(  1 30  )  /  IC'JRX(  130)  /  ISUMC/  ITSINV/LISTL/NBRANC/ 

1  KS/NP/JPCnL(ll )/ JPROW( 11 )/JPELEM(ll)/MAXC(10)/MAXC2(10) 

C 

1.A3TC  =  J*COl(NP*l)  -  1 
IF  (ITCNT  .Ct.  ITRFRQ)  GO  TO  20 
le  ( l F S A S  . u'«i.  0)  GO  TO  10 
RM TK  (21/1) 

1  format  (/•  OPTIMAL  INTEGER  SOLUTION  ') 

On  ‘.00  r  =  l/NP 

WRIT?  (21/2)  I 

2  FOP  ?AI  C  SUBPROBLEM'/ 14) 

J<FG  =  JFCOL(I) 

JMD  =  JFCOL(I  ♦  1)  -  1 

WKITP  (21/3)  (INCUMH(J)/  JsJPEG/JEND) 

3  FORMAT  (151b) 

♦or  conti  nu: 

4W£TL  ( 2 1  / *t )  I*CVAL 

s  FORMAT  ('  .i*X  OBJECTIVE  VALUE  ='/  16) 

PFTURN 

10  ^ITF  (2a/G)  INCVAL 

!»  F-WAT  (/*  Np  FEASIBLE  SOLUTION  FOUND  WITH  OBJECTIVE  VALUE  >'/ 
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1  110) 

:  "T’JK  N 

.(*  t.;> it:.  (2i/iCC)  itch? 

IOC  FQK'!AT(/*  SIMPLEX  ITERATIONS  =  */I8/':  COMPUTATIONS  TERMINATED.' 
IF  (IF i.AS  .KC.  0)  (id  TU  10 
WPITF.  (21/101) 

101  FORMAT  ('  bKST  INTEGER  SOLUTION  FOUND  IS:') 

UJ  600  l=l/NF 

WKITP  (21/2)  I 
JHFG  =  JFCOL(I) 

JF  ‘IL  =  JFCOLd  ♦  1)  -  1 

WRITE  (21/3)  (INCUMil(J)/  J=JBEG/JEND) 

COC  CONTI  n'JF 

a'BITF  (2i/o)  INC VAL 

<  FWATC  '<AX  n»JECTfVF  VALUE  DISCOVERED  =  '/I6) 

•  - TUr  1 
:•  NO 
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